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Use of Vectors in Vibration Measurement 


and Analysis’ 


CHARLES C. KENNEDY anp C. D. P. PANCUT 


Consolidated Vultee 


SUMMARY 


The authors’ main concern with vibration measurement and 
analysis is the identification and description of the normal modes 
of vibration of an airplane for the purpose of checking calculated 
modes to be used in flutter analyses or in the computation of 
dynamic loads. A technique for using four vibrators to separate 
normal modes directly while testing is presented. The basis of 
the technique lies in the fact that a single normal mode can be 
immediately recognized as such by its property of having the 
same phase response relative to the force at all points on the air- 
plane. The relation, between responses is either 0° or 180° and 
remains such if the vibrator frequency is changed. The procedure 
consists of systematically sweeping the point of application and 
frequency of the applied force until the responses are found to be 
exactly in or out of phase. The method is absolutely valid for 
separating two symmetrical and two antisymmetrical modes 
having resonant frequencies close together. 

For more complex systems, a technique is given for the identi- 
fication of normal modes on polar plots representing the variation 
of the response vector with forcing frequency. The basic element 
is that such a plot for a single mode is a circle; as the vibrator 
goes through a certain resonance, the tip of the total response 
vector describes an arc approaching that of the pure mode re- 
sponse. Identification of these arcs gives the normal modes, 
which might otherwise be confusingly located or completely hid- 
den in the current practice of considering just the amplitude of 
the forced response. 

Miscellaneous new points include treatments of the specific 
disadvantages of rotary vibrators, and the possibility of increas- 
ing the resonance sensitivity of the response measurements by 
differentiation with respect to the forcing frequency. Many 
other details of vibration analysis and record interpretation are 
given. 


INTRODUCTION 


og BASIC OBJECT of the paper is to make apparent 
the benefits of thinking in terms of time and space 
Received November 21, 1946. 
*The research on which this paper is based was done while 
both authors were Flutter Engineers at the Cornell Aeronautical 
Laboratory (formerly Curtiss-Wright Research Laboratory). 


t Now, Dynamics Group, San Diego Division. 


Aircraft Corporation 
U f 


vectors and to awaken a phase consciousness on the part 
of those who desire to bring about correlation between 
measured and calculated modes. The theory presented 
is a major step in the direction of such correlation; the 
main obstacle at present is the need for suitable vibra- 
tion exciting and measuring equipment to apply the 
theory. 
there is a demand, for elements satisfying the require- 


The designers can meet this need as soon as 
ments have already been developed in other fields. 


DESCRIPTION OF MOTION 


Normal Modes and Coordinates 


Since the main object of this paper will be that of 
determining and identifying the normal modes of vi- 
bration of mechanical systems, one must first consider 
what are normal modes and normal coordinates. 

It is a well-known theoretical conclusion that in the 
case of the forced oscillation of an undamped structure 
(or electrical network) the total motion is a sum of re- 
sponses in its characteristic modes of vibration; the lat- 
ter are variously referred to as the normal, principal, 
natural, or coupled modes of the system. The general 
problem of conservative systems has been studied by a 
(For 
example, see treatment by von Karman and Biot.') 
Rather than to attempt a complete and rigorous mathe- 
matical definition, only the pertinent properties of 


number of authors in various texts and papers. 


normal modes will be given here; they may serve as a 
definition for the nonmathematical reader. 

(1) Each normal mode responds to an applied force 
as a single-degree-of-freedom system—i.e., there is no 
coupling between normal modes. 

(2) In a normal mode, each point of the system os- 
cillates about its equilibrium position along a certain 


line in space, fixed relative to the equilibrium position 
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and straight when the oscillations are small enough so 
that all angles are equal to their sines. 

(3) In the case of simple harmonic vibration in a 
normal mode, all points move either exactly in or ex- 
actly out of phase with each other. That is, all points 
reach maximum departures from their equilibrium posi- 
tions at the same instants. 

(4) The shape of each normal mode is fixed for a given 
system and is independent of the magnitude, frequency, 
or location and direction in space of the applied external 
forces or of the deflections in other normal modes pres- 
ent. That is, in any given normal mode, the ratios of 
the deflections at all the points of a structure to the de- 
flection at an arbitrary reference point are always con- 
stant, and the directions of these deflections are fixed 
in space. These relative magnitudes and directions 
in space are characteristic of the normal mode, and their 
specification for every point in the structure will be 
referred to as the description of the ‘‘shape’’ of that 
mode. 

Actual measured airplane vibrations may rarely 
seem to have the properties of normal modes, but that 
is usually because of inability to separate normal 
modes when several are present simultaneously and not 
necessarily because the normal mode theory is inappli- 
cable. The properties of normal modes hold rigorously 
only for zero damping, which is, in general, closely 
enough approximated by aircraft structures. Damping 
is always present in actual systems and, as such, con- 
stitutes a controlling factor on the response of a normal 
mode as a whole; because the damping is small, how- 
ever, the mode shapes, phase relations between motions 
at various points, and couplings between the various 
degrees of freedom will be assumed to be unaffected. 

Turning now to coordinate systems, the following is 
quoted from von Karman and Biot:' ‘‘Let us consider 
a system of r degrees of freedom; for example, let us 
assume that the configuration of m mass points is given 
by r independent parameters. We call them the 7 
generalized coordinates qi, q2, .» gr and assume that 
we are able to express the Cartesian coordinate of the 
n points in terms of the 7 quantities qi, q, .» Yr by 
3n relations of the form: 


oe (q1, re qr) 
Ye = Vi (Qty Gor - ~~» Yr) (1) 


Zi (q1, q2; sey qr) 


i 
o 
| 


ll 


where 1 = 1, 2, ..., m. Then we are also able to ex- 
press the displacements 6x1, 61, 621; 6X2, dye, 420; 
vat ; 6xn, dyn, 6Z, in terms of the variations 6q:, 5q2, 


, 6g, of the generalized coordinates by differentia- 
For instance, we have: 


tion of Eqs. (1). 
5x; * Be (Ox ;/Ogx) 5qx ae 5X ix (2) 
= 1 = 


where 6x is the displacement at the point 7 in the co- 
ordinate g, or due to the mode k. We obtain analogous 
formulas for dy, and 62;. 


AERONAUTICAL 


SCIENCES—NOVEMBER, 1947 


A particular kind of generalized coordinate p, can 
be associated with every normal mode k; the more 
significant properties of this type of coordinate are as 
follows: 

(1) A normal coordinate deflection 6p, in the th 
normal mode determines the deflection (in Cartesian 
coordinates x, y, 2 or in terms of any generalized co- 
ordinates qi, gz, ..., g) corresponding to that mode for 
every point 7 in the system, the deflection relations 
being given by the mode “‘shape.”’ (Since the general- 
ized coordinates, and particularly the normal coordi- 
nates, are not commonly used to describe the location of 
the points but only their displacements from equilib- 
rium positions, the symbol 6 is often omitted below. ) 

(2) When using normal coordinates the total kinetic 
and potential energies of a conservative system can be 
expressed as sums of the energies in the various normal 
modes, where the kinetic and the potential energies of 
each normal mode can be calculated as the kinetic and 
the potential energies of a single mass.' There is no 
coupling between normal modes and, consequently, the 
equations of motion are greatly simplified, each equa- 
tion involving only one coordinate and its time deriva- 


sé 


tives. They are simply: 
Mp: + Kypr = F, (3) 
where k = 1,2, ...,7, Pr is the second time derivative 


of p,, and M,, K,, and F, are defined under (3) and (4) 
below. 

(3) In connection with every normal mode, a general- 
ized mass M, is determined so that 1/,p,? will equal the 
total kinetic energy of the structure in that mode—.e., 


M,p,? = > M(Xix? + Vin? + Six?) (4) 
t=1 
where 
d (3 ) Ox; 
Sn = —(— bp) = ba 
*  dt\ Op: Pr) = Pa op: 
d (2: ) ov; 
ae! me eo 
me ap, Prop, 
and 


ee oe d & ) —— Oz; 
BN\Oh dell lt Op: 


are the x, y, z components of the velocity of point 4 in 
mode k. A generalized spring stiffness K, is also de- 
termined by the relation wo, = VK,/M,. (w, is the 
natural frequency in rad. per sec. of mode .) 

(4) Corresponding to an actual force f that may be 
applied by a vibrator in a certain direction in space ata 
point 7 of the structure, a generalized force F; in the 
kth mode can be defined by the coordinate system if 
such a way that the work done by the force over a dis- 
placement at 7 in that mode, 5p, equals the product of 
the generalized force times the generalized normal dis- 


placement 6p,. Thus, 
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f(0q:/Opx)bpx = frbxx + fyivn + fbtu = Fidpe (5) 
where gq; is a linear coordinate measuring the displace- 
ment of point 7 along the line of action of f; f,, f,, f, are 
the components of the force in the x, y, z directions; 
and 
Ox; Oy; 02; 


i 
6 by 6 in = 6 ’ 6 | 6 
op Px Y ik ob. Px Six OD. Px 


When forces are applied at different points, 





OX iK = 


t 
Ox; Oy; 02; jo 
F, = Do (Sein + 2 fn + HS) (6) 
the summation being taken over all the points and di- 
rections in space 7 at which vibrators are applied to the 
structure. 

A generalized force may be made 0, thus eliminating 
the response in the corresponding mode by: (1) locat- 
ing a vibrator at a point of the structure which is an 
“absolute node”’ in that mode—i.e., 0x;/Op, = Oy;/ OP, 
= 02,/0p, = 0; (2) orienting a vibrator in space so 
that the force applied is normal to the response in that 
mode at the point where the vibrator is located—i.e., 
if 9 is the angle in space between the direction of the 
force applied f and the response in that mode 6p, at the 
point 7 where the vibrator is located 


cos §@ = (fr5X ix + f/Vx + S.52ix)/fop x = 0 


(3) locating multiple vibrators so that their effects 
cancel (e.g., use symmetrical locations to eliminate un- 
symmetrical modes); (4) varying the ratio of forces 
applied by the vibrators until the resultant is either 
located at an ‘‘absolute node’”’ as in (1) or at a node in 
the corresponding direction in space as in (2). (It 
should be pointed out that the preceding also applies 
without any essential modifications to the case 
when a vibrator is used which applies a torque and 
thus tends to produce angular rather than linear deflec- 
tions.) 

Since the generalized forces depend on responses, 
which are not the same for different modes, a single 
vibrator will apply different generalized forces in vari- 
ous modes. If two or more vibrators are not all in 
phase, the total generalized force F, can be found by 
vectorial addition of the component generalized forces 
from each vibrator with their respective phases, thus 
obtaining a new phase ye for each F, relative to a refer- 
ence vibrator. A rotary vibrator can be considered the 
equivalent of two linear-type vibrators contained in its 
plane of rotation, applying forces that are 90° out of 
phase along directions in space perpendicular to each 
other. Hence, a rotary vibrator will apply a general- 
ized force of a certain phase y, in the kth normal mode, 
but that phase will be different for different modes. At 
the instant when the vibrator mass is in line with the 
Projection of the normal mode deflection on the plane 
of rotation, the generalized force is a maximum, whence 
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Fic. 1. Phase yx of generalized force from rotary vibrator. 


the phase between any generalized force and the vibra- 
tor rotation can be determined (Fig. 1). 


Motion of a Point 


A displacement 6g, in a generalized coordinate k 
determines the displacement 6g, at the point 7 with 
components 6x, dix, 62 along the x, y, z axes (Fig. 
2). 6g is a vector quantity since it has magnitude and 
direction. It will be referred to as being a space vector 
since its direction is given in a space coordinate system. 

If 6g varies along the fixed direction sinusoidally 
with time, it may be represented as a different type of 
vector in a polar coordinate system, where the length of 
the vector is the amplitude of the motion and its direc- 
tion is given by the phase angle (wt + ¢). This type 
of vector will be called a time vector. The motion 
could be expressed analytically as 


dg = (5qu)4 Sin (wt + ¢) 
or 
du = (5qu)ae lies 
and the latter terminology will be used subsequently for 
convenience in mathematical operations, vectors being 
plotted in the complex plane. 

The x, y, z components of a sinusoidal motion along a 
straight line are also sinusoidal functions of time and 
have the same frequency and phase angle; the direc- 
tion cosines of the space vector will obviously remain 
constant throughout the cycle. It can now be seen 
that the condition of motion of a certain point along a 
straight line in space can only exist: (1) if there is but 
one normal mode present; (2) if the deflection at that 
point corresponding to each mode present has the same 
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Fic. 4. Calculation of h and a. 
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direction in space; or (3) if the responses in all the 
modes are in phase. 

When responses in two or more modes that do not 
have the same space and time vector directions are 
present at the same point of a structure, the motion 
along each normal coordinate direction can be resolved 
into its x, y, and z components of the form 

dx, = dx,,e! + 
foreach mode k. On any of the axes these components 
could be represented as time vectors and thus could be 
added together; the resultant time vector has a new 
phase angie, which will be different for the resultant on 
each axis. E.g.: 


r 
; + 
éx = >> bx,,€°% * Ov 
k=1 


—_ bx ei! + rz) 
~ he 


(7) 

When the x, y, and z components of the resultant 
have different phases, the motion cannot be a straight 
line in space because the direction cosines Ox/0q, 
Oy/Og, and 0z/Oq will vary with time. The motion 
is instead an ellipse, which is the most general type of 
motion for a point as long as the vibration in all the 
modes present is at the same frequency, as in the case 
of forced vibration without harmonics. 

The phase relation between the components of mo- 
tion along any two directions qg; and g2 can be deter- 
mined if the shape and orientation of the ellipse are 
known—a well-known technique with Lissajous figures. 
Consider the projection of the ellipse on the plane con- 
taining directions qi, g2 as in Fig. 3. The components 
of motion are: 


gi = dia COS (wt + 1) 
G2 = qra COS (wt + 2) (8) 


II 


The amplitudes of the two coordinates are qi, and qu 
the distances along the axis in question to the tangent 
parallel to the other axis. From Eqs. (8) it follows 
that gq; equals gi4 when wt = —q,. At that instant: 


(Q2)tor g = aa = 2a COS (—d1 + 2) (9) 


Hence, the phase relation between the two motions is: 


(d2 — ¢:) = + cos™! [ (G2) tor gi = qiA goa | (10) 


(Note the ambiguity of sign.) 

A linear vibration pickup measures the component 
of deflection along a certain direction in space which 
shall be called the pickup axis. When the point 7 to 
which it is attached oscillates sinusoidally along 4 
straight line, this motion is obviously in phase with its 
projection on the pickup axis. Corresponding to 4 
deflection p, in a normal mode k, the pickup measures 
dix = (0qi/Op,) px, Where q; is a linear coordinate meas- 
uring the displacement of point i in the direction of the 
pickup axis and 0q,/0, may be called the sensitivity of 
the pickup /ocation to mode k. The sensitivity of the 
pickup itself, S (= output reading /input), must include 


— 


> 


; 
i 





m 
sp 


mc 
cal 
for 


whi 
Ne 
thre 
ject 
mag 
whe 
am 
the f 
Le., 
the | 
xi 
two 
defle 
in Fj 
origit 
bzp re 
Hence 
(d2p , 


it is { 


ll the 


o not 
is are 
1otion 
solved 


nents 
ild be 
a new 
int on 


(7) 


ultant 
raight 
dx / Og, 
notion 
ype of 
all the 
1e case 


of mo- 
deter- 
yse are 
figures. 
1e con- 
onents 


(8) 
and qos 
angent 
follows 
ant: 


(9) 
ions is: 


| (10) 


yponent 
» which 
int 7 to 
along 4 
with its 
ig to 4 
1easures 
re meas- 
n of the 
tivity of 
y of the 
- include 








~-VIBRATION MEASUREMENT AND ANALYSIS 607 


phase as well as amplitude response; therefore, it is a 
time vector quantity. The overall mode sensitivity 
of the pickup, S(0q;/0p,), must be considered to be a 
vector quantity both in space and time. When the 
point of attachment moves in an ellipse, the projection 
of the motion on the pickup axis is sinusoidal, but it 
will have any phase angle depending on the relative 
orientation in space of g;—i.e., of the pickup axis—and 
the ellipse. 


Relations Between Motions at Various Points 


In a complicated structure such as an airplane, all 
parts participate in a normal mode vibration, and rela- 
tions between all the parts must be considered. When 
only one mode is present, some points will have no 
motion at all; they may be called absolute nodes. For 
any point one may speak of node directions in which 
there is no motion, and for any direction one may speak 
of node points at which there is no motion in that direc- 
tion. A linear pickup or vibrator located at a point 
and oriented in a direction in space in which there is no 
motion in a certain mode will not respond to or excite 
that mode. 

If more than one mode is present, the components of 
motion along fixed directions in space at various parts 
of the structure will have various phases, just as the 
x, y, components of motion at any one point will have 
different phase relations. Calculations of relative 
motions will thus have to be made vectorially.* Some 
special examples are now treated below. 

The method for the vectorial calculation of h and 
a, the bending and torsional components of the wing 
motion, from responses measured by two pickups lo- 
cated at the same spanwise station, makes use of the 
formula: 

6Zr — bZp 


a = tana = —— (rad.) 
Xp — Xp 


(11) 


which holds at any instant for small angles (Fig. 4A). 
If pz and dz vary sinusoidally, a does also and all 
three may be represented by time vectors. The pro- 
jection of the vectors on the real axis represents the 
magnitude of each quantity at the instant in the cycle 
when the vibrator force chosen for a phase reference is 
amaximum. If we combine the vectors according to 
the formula, the projections will also obey the formula— 
ie., if the formula holds for the projections, it holds for 
the vectors. 

Given the vectors for 5zy and 6zpg, the deflections at 
two points on the chord, the vectors for 6z, and a, the 
deflection at any point, and the wing twist are found as 
in Fig. 4B where éz, and 6Zp are laid off from the same 
origin O. The vector from ézy to dzp, (zg — dz) or 
ize relative to 5zy differs from «@ only by a scale factor. 
Hence, to get @ in radians, the magnitude of vector 
(62 — 62y) is measured to the scale of 5zr and dz, and 
itis then divided by (xg — xy), the distance between 
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the front and rear pickups. The formula for a may 
be rewritten a = (62, — dzy)/(x — xy) and solved for the 
deflection at any point, giving: 


dz, = (x — xp)a + dzp (12) 


Hence, once the vector for a has been obtained, a linear 
scale of x may be laid off on it and the deflection for 
any other point on the chord may be found, such as the 
trailing edge or a point on the aileron. Note that in- 
stead of a true node there is only a point of minimum 
amplitude at which the deflection 62,,;,. is always 90° 
out of phase with the angle of attack. When there is 
no fore-and-aft motion, the blur that the moving chord 
makes, as seen by an observer looking normal to the 
plane of motion (i.e., the projection of the envelope of 
the successive positions of the chord), is a hyperbola. 
Measurements made on this hyperbola may be used to 
obtain the phase relations between the motions at any 
points along the chord or between the angle of attack 
and the motion at any point. 

The angle of attack of the aileron may be found by a 
similar method, setting up the equation for linear 
measurement which will hold at any instant and apply- 
ing it to the vectors representing the sinusoidal varia- 
tion of the various deflections. 

It will often be useful to break down responses meas- 
ured at symmetrical points into components of symmet- 
rical and antisymmetrical motion. The normal modes 
of a closely symmetrical airplane are of two types. One 
type can be separated from the other because the fol- 
lowing mode shape data are known or can be assumed 
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in advance: (1) Ina given symmetrical mode, the ver- 
tical deflections at any symmetrical points R and LZ are 
equal (6zgs = 62,5); (2) in an antisymmetrical mode, 
these deflections are equal but opposite (6zg4 = —6zz4) 
The total deflections are: 


62, => O2Ls + 6224 


OZp = O2rs — b2ra 


(13) 


These two equations may be solved for the components 
in the two modes: 


62s = 1/5(62z + 5Zp) 


62r4 = 1/2(52z, — dzp) 


(14) 


The graphical construction is given in Fig. 5B. 
VECTOR RESPONSE TO FORCED VIBRATION 


Single Mode 


Eq. (8) above needs but the inclusion of damping 
effects to be applicable to actual systems. In airplane 
structural vibration work it is common practice to 
consider each mode k to have a structural damping 
coefficient g,, representing a damping force proportional 
to the spring force but 90° out of phase with it. The 
equation of motion in one mode is (dropping mode 
subscripts) : 


Mp + K(1 + ig)p = Fe (15) 


for a simple harmonic force.‘ (F, the amplitude of the 
force, may be complex.) Assuming that p = p,e™ 
is a solution and substituting, 


—Mw*p, + K(1 + 1g)pa = F 


Fig. 6 is a graphical representation of the relation in 
Eq. (16) [or in Eq. (15) if the vectors are considered to 
rotate with a circular frequency w]. 

From Eq. (16) the force required per unit displace- 
ment is obtained: 


F/p, = (K — Mw’) + iKg (17) 


and the variation with forcing frequency is represented 
vectorially in Fig. 7b for g = 0.25. The locus of the 
tip of F/p, is obviously a straight line parallel to the 
real axis, because the imaginary component of the 
force required for unit displacement is always equal to 
Kg. 

The plot of displacement per unit force amplitude 
(Fig. 7c) is a circle, with its ‘‘resonance diameter”’ 
(diameter through the resonant frequency point) per- 
pendicular to the generalized force. (This circle is the 
inverse conformal transformation of the straight line 
of Fig. 7b'.) Analytically, substituting the frequency 


ratio B = w/w,, where w, = VK/M is the natural fre- 
quency, one gets: 


Pa 1 _ 1 
F  (K — Mo") + iKg K[(1 — 6?) + ig] 


(16) 





(18) 
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Nondimensionalizing the displacement in terms of the 
static deflection F/K, the following expressions result: 


Kpas/F = 1/((1 — 8?) + ig] (19a) 





or 
Pa _ = : §tan~! [—g/(1 — 8?)] 
F/K V(1 — 63? + g?- (19b) 
Kp,/F = ue’ (19¢) 


in polar form, and 


1 — 8 





a a a 
F/K LO — 6)?+ 2 +4 aG-#) +e (19d) 


Kpa/F = ur + tpy (19e) 


in rectangular coordinates. (The symbols introduced 
are defined by comparison of terms.) 

The response at a point 7 and along a direction in 
space or coordinate q;, induced by an exciting force of 
the form F cos wt is: 


dq: F , 
ig, = = - pees Gt + 6) (20a) 


DK 
where the phase angle ¢ is the lead of the motion rela- 
tive to the generalized force. The complex vector is 
represented in a coordinate system rotating at a fre- 

quency w as: 

sq, = 5 = uei® (20b) 
The generalized force corresponding to a single linear 
vibrator applying a force f at a point 7 in a direction 
gq, is F = (0q,/0p)f by Eq. (5). The response induced 
by this force along direction q, at a point 7 is obtained 

by replacing in Eq. (20b): 
- 09: 0g, f ip 91) 
ace at il " 
It can be seen that the response at 7 due to the vibrator 
force at j is equal in magnitude and phase to the re- 
sponse that would be induced at 7 if the same vibrator 
were located at i, because the same partial derivatives 
would be involved (Maxwell’s reciprocity theorem). 
(In the preceding statement i and j should be under- 
stood to represent points as well as directions in space.) 
The polar form of Eq. (19c) is familiar, but the rec- 
tangular form with the new symbols for the real and 
imaginary compoents yr and yu; of the nondimensional 
response is unusual.’ The usefulness of these terms 
will be seen later. Figs. 7d and 7e give plots of », 4, 


br, and yu; vs. 6? for a damping coefficient g = 0.03. | 
It should be noted that resonance is indicated by max: | 


ima in pw and yu; (their derivatives with respect to 6? are 
0), and by points of inflection in ¢ and ye (maximum 
slope or maximum derivative with respect to 6°). 
Also pp = 0 and ¢ = 90° at resonance. 
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Fic. 7. Single degree of freedom response, 


For a rotary or inertia vibrator with a constant un- 
balance, the force does not have a constant amplitude 
since the forcing frequency is varied; it is equal to the 
unbalance mass m,, multiplied by the absolute accelera- 
tion of its center of gravity—i.e.: 


Actual force = —m,(gr + Gr) = m,(Rw*e"* — go) 


where jg = acceleration of m, relative to v, the point 
of attachment to the structure; g, = absolute accelera- 
tion of the point of attachment; R = arm of rotation; 
and w = angular frequency. The equation of motion 
in any normal mode is of the form: 

Mp + K(1 + ig)p = (dq,/0p)m,(Rw°e™ — G,) (22a) 


which, after rearranging and substituting i = 
(0g,/p)p, becomes: 


lp)? ‘ ; 

M + m,(0q,/0p) b+ a = m,Rete (22b) 
Og,/Op 0q,/Op 

This represents the new system shown in Fig. 8a for 
which: the generalized mass is increased by the con- 
tribution of the vibrator mass; the generalized mass 
and spring are changed in proportion to the normal 
mode deflection at the point of attachment and the 
driving force is directly proportional to the frequency 
squared. 





Again assuming a solution of the form p = pye™, 
the response is obtained as: 
Pa -_ m,(0q,/0p)R uBe'® (23) 


~ M+ m,(0q,/ Op)? 


where 
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Fic. 8. Frequency response plots for various dampings. 


p= 1/V(1— B2)2?+ 8% 9 =tan—[—g/(1 — B)] 


and 
B = w/VK/[M + m,(2p./dp)*] 


For the case of a rotary vibrator, an additional phase 
angle y will be added in each mode as explained with 
Fig. 1. 

Figs. 8b and 8c are vector plots of (1/u8?)e 
the inverse, u8%e'*, for various values of g. These may 
be considered as nondimensionalized graphs of, respec- 
tively, unbalance required per unit displacement and 
displacement response to a unit unbalance for variable 
forcing frequency for the system of Fig. 8a. They may 
be also considered to be plots of the force required for 
unit acceleration and the acceleration response to a 
unit force for the system shown in Fig. 7a, since: 


~*? and 


. 
F/M 


paw” 


paw” 9 1 
wa" nae Oe OO 


 Feog?/K 


Fig. 8b is also, by coincidence, a straight line, sometimes 
called the ‘‘damping line”’ in graphical flutter analyses.°* 
Similar curves corresponding to products of ye’® with 
any other powers of 8 are more complex. (For systems 
with viscous damping, only the plot of force required 
per unit velocity is a straight line;® the curve of force 
per unit displacement is a parabola.) When damping 
is-small, however, a nearly straight line is always ob- 
tained in the vicinity of the natural frequency. Note 
the simplicity of Figs. 8b and 8c as compared to the 
usual separate phase and amplitude response plots. 
The effect of lowering the damping, which is that of in- 
creasing the diameters of the circles and of expanding 
the frequency scale, can be readily observed. For 
any one given damping, the frequency scale is most 
expanded in the vicinity of resonance; as shown be- 
low, this property is the basis of a novel method of lo- 
cating natural frequencies of considerable utility when 
the curves are complicated by the presence of many 
modes. 
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Vector Derivatives of Response Vectors 


Before considering the rate of change of the response 
curves with frequency ratio, a brief review of the differ- 
entiation of vector quantities in general will be made. 
In Fig. 9a let curve C be the path of motion of a point 
P constrained to move in the XOY plane. The posi- 
tion of the point is given by the tip of the vector V(t), 
which is a function of a certain parameter ¢. The de- 
rivative with respect to ¢ is: 


aa dV. . | [V(t+ At) — vO] i 
, =o = —_—_—_———— a ) 
V(t) r limit | a (25) 


At — 0 


(The minus sign indicates vector subtraction.) V’(t) 
is a vector itself, having the direction of the tangent to 
curve C and a magnitude equal to the “‘speed’”’ of the 
tip of V(t), ds/dt, along the path C. The curve C’ is 
the hodograph of C. 

Applying this to differentiating ye with respect to 
8?, it may be shown that 


(ue’*)’ = (ue'*)? = pre? (26) 
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i.e., that the square of yue'® equals its first derivative 
with respect to 6. This can be proved by using the 
properties of the circular vector plot, Fig. 9b. The 
angle of the tangent is 29, which is the angle of the de- 
rivative. The length s of the arc on the curve is the 
radius times the angle subtended; hence 


ds = R X 2do = dd/g 
The argument of the derivative is: 


ds 1 do 1 d E _ ( —g )| . 

= — Se a — = toad 

d(g*) g d(B’) g d(?) 1 — 6 4 
Q.E.D. 





Fig. 10a gives the locus of the tip of u?e"?* in the com- 
plex plane. The equation for this curve in polar coordi- 
nates, obtained by eliminating 6? between yu? and 2¢, is 


mu? = (1/2g")[1 — cos (2¢)] (27) 


On plots of the amplitude and phase (u? and 2¢, re- 
spectively) of the derivative vs. 8”, resonance would 
again be indicated by a maximum or an inflection point, 
respectively, in the two curves. The peak in the yz? 
plot, however, is sharper than that in uw, and the maxi- 
mum slope of the 2¢ curve is greater than that of ¢; 
consequently, it may be said that the derivative is 
more sensitive to a resonant frequency than is the 
original quantity. 

In theory, the vector differentiation with respect to 
8? of ye’* can be repeated any number of times and the 
nth vector derivative can be shown by inductive reason- 
ing to be: 





e" ; ' . 
Haye We) = at wh Piel +O (28) 


Each successive differentiation will sharpen the peak 
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[d/d(A*)] (ue?) = (uete)’. 
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in the curve of amplitude response vs. 6? and also 
steepen the inflection slope of the curve of phase re- 
sponse vs. 6° near the natural frequency; in other 
words, the sensitization to resonance is increased by 
differentiation. 

The same effect can be seen in another manifesta- 
tion. The real and imaginary components of the deri- 
vative are equal to the derivatives of, respectively, 


Response for two modes, natural frequencies nearly equal. 


the real and imaginary components of the response 
vector. In Fig. 10b are given plots of up’ and u;’, the 
rates of change with §? of the real and imaginary com- 
ponents of ye’*: 
d 


tr’ = d(8?) (u cos >) _ we cos (2¢) 


d 
(ur) = is’) 


and similarly 
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My’ = uw’ sin (2¢) 


These derivatives too have greater slopes and sharper 
peaks than the original quantities. 

The importance of the derivatives in practical vibra- 
tion work lies mainly in giving another point from 
which to view the vector plots. However, more than 
one differentiation will be impractical with the present 
uncertainty of vibration measurements made on me- 
chanical systems. It might find more direct applica- 
tion in electrical problems where experimental inac- 
curacies may be smaller. 


Two Modes 

To find the response 6g; at a point 7 and in a pickup 
direction g; when two modes are excited simultane- 
ously, Eq. (2) can be used: 


Og: Og: 
6g: = =— 6f1 + — bp 
dt dp: 1 Dp, 2 


If the oscillation is harmonic, Eq. (20b) applies for each 
mode: 


(29a) 


6q; = que’* = ous fs 


Og: Fo 
) D1 K, Mi Me 


t(gi + ya) i: 
+. Op. Ke 


ei (* + y2) 


(29b) 


where y is the phase lead of the generalized force rela- 
tive to some reference such as the phase of the motion 
of the vibrators. _ 

A typical example of such a summation for two ar- 
bitrary modes is shown in Fig. 11, which is a polar plot 
of the locus in the complex plane of the tip of vector 
qe’*' as given by Eq. (29b) above, when the forcing 
frequency is varied from 0 to ~. As can be seen, it is 
a curve having two nearly circular parts corresponding 
to the two modes. This vector plot could be obtained 
in three different ways: (1) adding separately the real 
and imaginary components of the responses in the two 
modes (from Figs. 1ld and 1le); (2) adding vectorially 
the responses in the two modes taking into account the 
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respective magnitudes and time phases (from Figs. 
1lb and llc); and (3) geometrically by finding the 
locus of the point midway between points corresponding 
to the same frequency on the two single mode vector 
plots. (This construction, exemplified in Fig. lla for 
N = 1475 cycles per min., gives the vector 0.5(g: + q), 
and, therefore, the single mode response circles have to 
be drawn with a double radius.) Note that the peaks 
and inflection points corresponding to resonance on the 
de,qi and the differentiated curves occur more closely 
to the natural frequencies of the normal modes than 
they do on the resultant g and ¢ curves. 

Figs. 12 to 16 are similar plots for other examples of 
responses with more than one mode present. Table 1 
lists the values of the various constants in Eq. (29b) for 
these figures. 

Fig. 12 shows the effect of having two modes with 
natural frequencies close together. The vector plot 
still shows two distinct circular parts, but the connect- 
ing loop in Fig. 11 has become merely a bulge; two 
points of inflection in the curve replace the loop. 

In Fig. 13 one mode has an extremely low damping, 
g = 0.01, which increases its relative amplitude. The 
smaller damping mode tends to obscure the mode with 
the larger damping, but the two circular parts are still 
apparent. 

Fig. 14 shows the variation of response with pickup 
location along the wing chord when two modes are 
present. The curves for points A and D exemplify 
the case when 0q;/0p, has opposite signs for the two 
modes which causes part of the curve to fall in the 
positive imaginary region.* For curve D the two 


* Only when a single linear vibrator or its equivalent is em- 
ployed, and then only at the point of application and in the di- 
rection in space of the exciting force, must the response vector 
plot be contained entirely in the negative imaginary half of the 
complex plane. This is due to the fact that all mode response 
circles are on the same side of the force phase reference test which 
is a corollary of the Driving Point Impedance Theorem.’ When 
single mode response circles fall on both sides of the real axis, 
the paradoxical situation may arise where a resonant frequency is 
indicated by a minimum rather than a peak in the amplitude 
response curve (see Fig. 16). 





TABLE 1 
Summary of Mode Data for Figs. 11 to 16 
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Natural ing a 
Fre- Coef- ogi Fr 
Fig. Mode quency ficient ope Kk vb 
11 1 1,400 0.03 2 0 
2 1,500 0.05 3 0 
12 1 1,500 0.05 3 0 
2 1,550 0.03 2 0 
13 1 1,400 0.01 2 0 
2 1,500 0.05 3 0 
14 — Point 
AO: B C O D 
1 1,400 0.03 ~ § 8.86 2 0.5 0 —1 : 
2 1,500 0.05 -3 0 36 7 9 ; 
15 1 1,400 0.03 2 ° 
2 1,500 0.05 3 
16 1 1,400 0.03 _ : 
2 1,450 0.03 —0.75 0 
3 1,500 0.05 3 ne 
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nearly circular parts are not even separated by two 
inflection points; they are only distinguishable by the 
increased rate of change with frequency of the response 
vector at resonance. The frequency scale is marked in 
constant intervals of 10 cycles per min. (as in all the 
other vector plots) to make the changes in ds/dN and 
correspondingly in ds/d(8*) apparent. In Figs. l4c to 
14h the various types of response for points A and D 
are compared. Note that only one peak in amplitude 
occurs in either curve (Fig. l4c), but the amplitude 
of the derivative shows peaks for both modes (Fig. 
14g). 

The vector plots of Fig. 14a are obtained as follows: 
If the responses in the two modes at O; and Os, the 
nodes of the modes for motion in vertical direction, are 
given, the respective circles can easily be constructed by 
Fig. 7. Then, at any frequency, the response vectors 
for these two points are combined, according to Eqs. 
(11) and (12) and Fig. 4B, to find the response at any 
other points A, B, C, .... along the chord. To illus- 
trate this, the ‘‘chord’’ lines (loci of the tips of the 
response vectors for various points along the wing 
chord) are drawn in for 1,400 and 1,500 cycles per min. 
If the construction is repeated for a number of fre- 
quencies, curves A, B,C, .... of Fig. 14a corresponding 
to points A, B, C, ... of Fig. 14b are obtained. (A 
different interpretation can be given by considering the 
curves O;, Oo, A, B, C, ... to represent responses at the 
same point P but along different directions in space. 


Response for two modes, generalized forces 90° out of phase. 


Particularly then O, and Oy» are the node directions— 
i.e., the directions of zero response in, respectively, 
modes 1 and 2 at point P.) The method can be used 
to find the response curves at any point in the structure 
and for any direction in space, because the points along 
the chord line (extended to infinity) represent all pos- 
sible combinations of proportions of responses in the 
two modes. That is, all the possible shapes of fre- 
quency response curves will be given by this construc- 
tion, although the actual magnitudes might be different 
at other parts of the structure. Note that by the reci- 
procity theorem the same figure would also represent 
the responses at the present vibrator location if the 
same vibrator were placed at points A, B, C, . .. . of the 
structure. 


In Fig. 15, the generalized forces are assumed to be 
90° out of phase to exemplify the effect of using a 10 
tary vibrator or, in general, linear vibrators that art 
out of phase with one another. As the resonance di 
ameters of the various response circles have to be per | 
pendicular to the generalized forces in the respectivé 
modes, the maxima in the imaginary components wil | 
indicate natural frequencies only when these general | 
ized forces are all in phase. In a case as that of Fig. 5 | 
resonances must be identified by the fact that they cor | 
respond to points of maximum ‘‘velocity’’ along th 
curve as N changes (i.e., points at which ds/d(s*) 84f 
maximum), 
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Fic. 16. Response for three modes, one hidden. 


More than Two Modes 


It is apparent that the vector responses at a point 
can be synthesized for any number of modes when 
given g, and (0g;/Op;,)-(F,/K,). The circles for re- 
sponse in any one mode may be displaced farther and 
farther from the origin, but they will usually still be 
evident, having distortions similar to those seen in the 
examples given for two modes. The phase and ampli- 
tude curves will flatten out as the circles are displaced 
farther, and, consequently, it becomes more and more 
necessary to make vector plots to identify modes. 

In some cases, however, the circles may be hidden. 
This may happen when the resonance diameter of the 
circle lies on the side of the reference opposite to that 


shown in most of the examples. In this case the circle 
may only accentuate a loop or inward bulge of the 
curve by increasing its size, and thus, adding no new 
characteristic to the plot, the corresponding mode 
might go unnoticed. An example is given in Fig. 16. 
Mode Nos. | and 3 are the same as shown in Fig. 11. 
A mode No. 2 has been added at NV = 1,450. Compar- 
ing these figures, the most apparent difference is the 
increase in the size of the loop; an unwary person might 
think the response represented only two modes. On 
closer inspection a criterion appears: comparing As for 
a constant AN three maxima are found in Fig. 16a as 
compared to only two on Fig. lla. This is obvious in 
the plot of the amplitude of the derivative, Fig. l6e, 
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but in none of the other curves is the 1,450 mode appar- 
ent. Again, the sensitivity to resonance of the deriva- 
tive is illustrated. (The previous derivatives were with 
respect to 6’ or, in effect, to N*. The shape of the deriv- 
ative with respect to N will be little different, since 
dq/dN = 2N[dq/d(N?)].) In extreme caseS, one dif- 
ferentiation may not bring out the hidden peak in the 
valley between two obvious peaks, but in theory it 
would eventually appear in the amplitude curve as 
higher derivatives were taken. In practice, however, 
the uncertainties in measurements are magnified with 
differentiation; hence, with the present inaccuracy of 
experimental data more than one such differentiation 
may be impractical. 

The following are the more important conclusions re- 
garding frequency response vector plot shapes: 

(1) In the vicinity of a resonance the shape ap- 
proaches that of a circle. 

(2) In most cases, for frequencies intermediate be- 
tween two resonances the response vector tip describes 
a characteristic loop such as that of Figs. lla, 13, 14a 
(points B, C), and 15a. 

(3) When resonant frequencies are close to one 
another (Fig. 12a), or when dampings are large or re- 
sponses are different in size and angular positions [Fig. 
14a (point A)]—i.e., whenever the effect in the vicinity 
of one natural frequency due to another mode is sub- 
stantial, no such connecting loop can be observed. 
However, the radius of curvature, which usually is 
directed toward the origin in the vicinity of resonances, 
will still be directed the opposite way in this transition 
region, the change being accompanied by two inflection 
points which replace the loop. 

(4) At resonance ds/d(6*) is always a maximum. 
This can be seen either on plots of this quantity (the 
amplitude of the derivative) vs. 6”, or then directly on 
the vector plots when equal intervals of Af or AN are 
laid off. 

All of the discussion in these sections can be applied 
to plots of acceleration response to constant vibrator 
force amplitude and closely to other forms of response 
and other types of damping, even in the cases when the 
latter is sufficiently large so that the vibration modes 
do not have the simple characteristics listed at the be- 
ginning of this paper. The departures of the various 
kinds of response plots from the circles discussed here 
will be negligible in actual airplane vibration practice. 

A similar type of vector representation of frequency 
responses is used commonly in electrical problems such 
as the oscillatory stability of feedback amplifiers.'!» 
More recently, the method of plotting a forced oscilla- 
tion response in the complex plane has been mentioned 
in connection with the design of servomechanisms.” 
The problem of airplane vibration is made more compli- 
cated than these by the fact that an airplane is a dis- 
tributed system with fixed and unknown stiffness, 
damping, and inertia constants, and the experimental 
technique and apparatus is less well developed.  For- 
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tunately, the damping of structures is usually low 
enough to enable approximations by normal modes, 
which is not always possible in connection with other 
kinds of systems. 


SEPARATING MODES 


The Problem 


We are now ready to discuss the main problem of ex- 
citing, measuring, and analyzing the vibration modes of 
a structure and of airplanes in particular. The object 
is first to “identify’’ the various normal modes in a 
certain frequency range and then to ‘‘describe’’ each 
mode of interest. 

The identification of a mode will be made by giving 
its natural frequency and if necessary, some shape 
characteristic that will distinguish it from other modes 
with resonant frequencies in that vicinity. Rough de- 
scriptions such as bending, torsion, wing, tail, or fuse- 
lage modes may thus serve a useful purpose; however, 
one must be careful not to assume that the mode shape 
is known when such a label is attached, because several 
of these uncoupled degrees of freedom may participate 
in an actual normal mode. It is safer to use just the 
natural frequency as an identification whenever pos- 
sible. : 

The description, to be complete, must include (1) the 
natural frequency, (2) the damping, (3) the generalized 
mass or generalized spring stiffness (one will determine 
the other), and (4) the complete shape of the mode. 
The shape will give the motion of every point corre- 
sponding to a unit displacement in that normal mode— 
i.e., it will determine Ox,/O0p,;, Oy,/Op,;, and Oz,/Op, at 
all the points of the system. [See Eq. (1).] 

To carry out the identification, the total forced vi- 
bration of enough points of the structure must be meas- 
ured so as to be sure that all modes will show their 
resonance characteristics in at least one of these re- 
sponses. For this purpose it is now common practice 
to use plots of the total amplitude response vs. forcing 
frequency. One of the main objects of this paper is to 
show that the best way—and, when the total response is 
especially complicated, the only way—is to plot the 
response vector (i.e., to take into account simultane- 
ously both amplitude and phase response). 

The first step in trying to simplify a complex total 
response is to eliminate as many undesired modes as 
possible by using vibrators that do not excite those 
modes. In general, with one linear vibrator it is pos- 
sible to eliminate at most three modes; the point of 
application of the force will have to be an absolute node 
of the first mode, and the direction of the force must be 
perpendicular to the plane containing the space direc- 
tions of the responses at that point in the second and 
third modes. In practice this is not always possible be- 
cause absolute nodes may be inaccessible. On the 
other hand, several modes may have common node 
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points or directions and then more than three modes 
may be eliminated (e.g., a vibrator acting in the plane 
of symmetry of a symmetrical structure will excite 
none of the antisymmetrical modes). A fourth mode 
may also be eliminated by adding a second vibrator so 
that it also will not excite the first three modes and then 
varying its force output until its contribution to the 
generalized force in the fourth mode cancels that of the 
first vibrator, Eq. (6). Similarly, each additional vibra- 
tor allows one more mode to be eliminated. 

The above paragraph applies identically to the elimi- 
nation of responses in undesirable modes by properly 
locating and orienting linear vibration measuring pick- 
ups. Forced responses at different points of a structure 
may be combined mechanically, and the pickup re- 
sponses may be combined electrically before recording. 
Finally, the recorded responses may be combined 
graphically by the methods discussed in connection 
with Figs. 4, 5, and 14. (It should be noted that re- 
sponses in undesirable modes are eliminated from all 
the points of the structure when vibrators are properly 
used, while only one response at a time may be thus 
“purified’’ by combining pickups.) 

For the identification of modes of a symmetrical 
structure such as an airplane, the most practical setup 
is that in which two or four linear vibrators are em- 
ployed. A technique for the separation of normal 
modes which can be used directly during the vibration 
test will be presented. 


Use of Four Vibrators (on Symmetrical Structures) 


The equipment must consist of: 

(1) Four linear vibrators that can be kept closely 
in phase under varying load conditions and may be 
connected so that two in a pair will have exactly equal 
or opposite forces. These forces must be variable but 
need not be entirely controllable (for example, they 
may depend to some extent on the response of the struc- 
ture); an accurate indication of both their magnitudes 
and time phases must be available, however. An abso- 
lutely immediate force ‘cutoff is extremely useful for 
taking decay records. Electrical “loudspeaker’’-type 
vibrators with heavy field coils acting as seismic ele- 
ments may be made to satisfy these requirements. The 
“voice” or “armature’’ coils of any pair of vibrators 
can be connected in series and the direct current fields 
are adjusted to equalize the forces applied. The alter- 
nating voltage across the voice coils is supplied by a 
beat-frequency oscillator, and the magnitude and phase 
of the applied force are obtained from a record of the 
armature current to which the former may be assumed 
to be directly proportional. (Such a system, originally 
developed by the Lockheed Aircraft Corporation, has 
been used by the writers with good results.) 

(2) A number of linear vibration measuring pick- 
ups (no less than four) with which the phase and ampli- 
tude responses at various points can be compared simul- 
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Fic. 17. Resultant force—two vibrators. 


taneously and directly with the vibrator force while 
testing. The phase sensitivities of the pickups must 
be identical, known, and preferably equal to zero. 
Ability to record when desired is almost essential for 
complicated systems and for determining the damping. 
At present, these requirements may be met more or 
less satisfactorily with electrical accelerometers and a 
recording oscillograph with a viewing screen. (The 
accurate measurement of phase responses on sinusoidal 
records is somewhat tedious, however.) 


The procedure is as follows: 


(1) Attach the vibrators at symmetrical locations 
on the airplane in.such a way that when the force out- 
puts of the right and left vibrators are either equal or 
exactly opposite; only symmetrical or, respectively, 
only antisymmetrical modes are excited. The system 
now reduces to the general case of two vibrators, and 
this will be assumed in the following discussion. To be 
more specific, let us also suppose that these two vibra- 
tors and two linear pickups are oriented in vertical di- 
rections and are located at the same spanwise station 
of a wing—e.g., at the front and rear spar on the same 
rib. (Possible alternative arrangements would be 
those in which either the two vibrators or the two pick- 
ups or both are located at two different points along the 
same wing spar and have parallel directions in space or 
then are placed at the same point of the wing but in 
directions in space at right angles to each other as, for 
example, in a chordwise and a direction normal to the 
plan form of the wing.) Locate the other pickups at 
any points and in any directions such as to best cover 
simultaneously all the modes of interest. At least one 
pickup should be used on the opposite wing to check the 
symmetry of the responses. 

(2) When two vibrators are located on the same 
chord, it is possible to change the point of application 
of the resultant force by varying the ratio of the two 
forces. With symbols and sign conventions defined 
as shown in Fig. 17, the position of the center of force 
is given by 


(30) 


C.F. = (xe — xr)/[1 + (F/R)] 
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where F/R is the ratio ot the front and rear forces. 
The generalized force in any mode will depend on F/R 
since it is proportional to Ogc. r,/Op;, the displacement 
at the center of force for a unit displacement in the 
normal coordinate. If the resultant acts at the node 
of a certain mode k—i.e., Ogc. p./OP, = O—that mode 
will be eliminated. 

If only two modes are present (symmetrical or anti- 
symmetrical), start by running the vibrators at the 
lowest frequency in the range of interest and vary F/R 
until the various pickup responses are observed to be in 
phase—i.e., until the resultant force is applied at the 
node of one of the modes. 

(3) Increase the forcing frequency until a peak am- 
plitude per unit force isfound. This is the first resonant 
frequency, and the corresponding normal mode may be 
identified directly and its shape measured. Damping 
may be obtained from a record of the response during a 
frequency sweep test or from a curve of free vibration 
decay from this forced condition. 

(4) With the vibrators running at the first resonant 
frequency, again vary F/R until a new in-phase re- 
sponse is found—i.e., put the center of force on the 
node in the first mode. Thus, the second mode will be 
excited alone. 

(5) Vary the forcing frequency until the peak am- 
plitude indicating the second resonant frequency is 
found, and measure as before. 

With the two following exceptions the procedure is 
absolutely valid for two modes. If the two modes 
have the same natural frequency and damping, all 
points would have the same phase response at all fre- 
quencies, and the node positions could not be found as 
above. Instead, the distinction must be based on the 
fact that the resultant mode shape would vary with 
the position of the resultant force. If the nodes of ver- 
tical motion in the two modes happen to coincide, it 
will not be possible to excite only one mode at one time. 
Another spanwise location could be tried. When many 
modes exist, the procedure is still closely valid either for 
separating two modes having natural frequencies close 
to one another but well-separated from other resonances 
or, in general, to eliminate from the total reponse 
the component corresponding to any one given 
mode. 

Decay records serve to check the purity of the modes 
thus found, giving an independent means of measuring 
the resonant frequency, mode shape, and damping. 
The response in each normal mode decays exponentially 
with a frequency that is the natural frequency of that 
mode. If responses in two or more modes are measured 
by any pickup, that will be manifest in the number of 
different frequency components in the recorded wave 
form. If only two modes with frequencies close to- 


gether are present, the record will be a beat that can be 
analyzed to give: the two frequencies and ampli- 
tudes. "4 
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Analysis of Records 


As the interference of the numerous modes gets 
larger with increased complexity of the system, re. 
course must ultimately be made to vector plots of 
response. The procedure is to: (1) record responses 
during preliminary frequency sweep tests with arbitrary 
vibrator locations (such as on nodes of calculated, or es- 
timated coupled, or uncoupled modes); (2) identify 
the modes from record analysis; (3) obtain complete 
mode descriptions from final forced and free vibration 
tests. In the last, use vibrators and pickups at locations 
chosen on the basis of the findings under (2), so as to 
best excite and measure the various modes of interest. 

Frequency sweeps can be recorded for one center of 
force or vibrator location, and then curves for various 
points along the wing chord may be found graphically 
from any two pickups by the method of Fig. 4B (as in 
Fig. 14). Again the reciprocity theorem is useful: 
These curves give the response at the position used for 
the vibrator if the point of application of the force were 
moved along the chord, thus making it unnecessary in 
many cases to actually move the center of force. 

The general method of identifying the modes from the 
frequency response vector plots consists of drawing the 
best circle possible through each circular arc on which 
maximum in ds/d(6*) (Fig. 9b) is distinguishable. This 
circle is the best approximation obtainable from the 
given data to the response measured by the respective 
pickup in the corresponding mode. The reason for this 
is that, when the vibrator is sweeping through a certain 
resonant frequency the tip of the total response vector 
will tend to partly describe the circle corresponding to 
that mode as the components due to other modes re- 
main more or less constant. The resonance diameter is 
drawn perpendicularly to the generalized force in the 
mode (i.e., to the force phase reference when in-phase 
linear vibrators are used), and it should pass through the 
point of maximum ds/d8 on the curve; at its other 
end, it will determine the ‘displaced origin’’ to be used 
when measurements are made. 

The record analysis is exemplified in Figs. 18, 19, 
and 20. Vector response plots from previous figures are 
analyzed as if they were actual record measurements. 
On these, the best single mode response circles are 
drawn through the “‘measured’’ points; subsequently, 
the natural frequency is found at the end of the reso- 
nance diameter, and a plot of amplitude (measured from 
the displaced origin) vs. frequency is made to find the 
peak amplitude and the damping. To obtain the latter, 
use was made of the relation that for a single mode with 
small damping: 





2q/qmaz. |N 
= —e  - 31 
Oe Vt = (g/Qmaz.) | Nn oe 

which, for ¢/qmaz, = 0.707, becomes: 
Pag 2| Bo.707 vet 1| (32) 
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where (707 is the value of the frequency ratio N/N, at method is obvious. The improvement in frequency 
the point where the amplitude q is 0.707 of the peak determination is not apparent from these tables; it 
amplitude gmar. or where the response vector is 45° should be noted, however, that in a case such as that 
away from the resonance diameter. in Fig. 19 in which no peak occurs, the ordinary method 
In the table given in each figure the results are com- would give no answer at all, unless exceptional benefit 
pared with those obtained by the usual method of of the doubt be given the analyst by assuming that he 
neglecting completely the effect of other modes present. would recognize the existence of a mode near the point 
The per-cent errors, based on the known true response of maximum curvature of the amplitude response. An- 
in the pure mode, are given for the two methods; nega- other such example is the 1,450 mode in Fig. 16, which 
tive errors indicate that the ‘‘measured’’ values are would be completely hidden in the ordinary analysis. 
too small. With regard to damping and peak ampli- Before drawing in the response circles, the construc- 
tion indicated in Fig. 14a may first be used to find the 
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vector plots for various points along the chord; the 
circles then may be fitted to the simpler of these curves. 
On a diagram such as that of Fig. 14a, note that, since 
the diameter of the circle response in a given mode re- 
duces to 0 at the node, the ‘“‘chord” lines will tend to 
pass through a common point for a frequency range 
near the respective resonance. The relative position 
of this point and that of the tips of the various response 
vectors along the ‘‘chord”’ lines will give the actual node 
location of that wing station. 


CONCLUSIONS AND RECOMMENDATIONS 


The authors have presented a theory for the measure- 
ment and interpretation of forced vibration test data, 


giving several new facts, analyses, techniques, and 
points of view. 

Since no experimental examples are given, one 
might well inquire as to the validity and use- 
fulness of the theory in actual practice. The experi- 
mental work done by the authors has included a thor- 
ough survey of the vibration modes of a single-engine 
airplane using single mechanical linear and rotary vibra- 
tors and also multiple electric exciters. Special atten- 
tion was paid to the wing vibration characteristics in 
the frequency range from 1,200 to 2,200 cycles per 
min. which presents a confusing picture on many air- 
planes of this type; this is due to the fact that there 
are between five and ten modes involved corresponding 
to wing torsion, second bending, and fore-and-aft 
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motion, aileron and flap control systems, landing gear, 
fuselage, and tail degrees of freedom.* 

The usefulness of the four vibrator technique de- 
scribed in this paper was well demonstrated by the 
fact that some pure decay records resulted, showing 
that the responses in all the modes but one could thus 
be effectively eliminated. A large number of frequency 
response vector plots were constructed, and mode 
circles were identified after graphical methods were 
used to simplify these curves. Despite relatively large 
uncertainties due to record reading, inconsistent per- 
formance of the exciting and measuring equipment 
and variation of the airplane condition over the pro- 
tracted time of testing, fair correlation between the 
various tests was possible. This analysis served the 
important function of proving the necessity for better 
methods of ground vibration and pointed out direc- 
tions for improvements. 

To get reliable data for such analysis, vibrators 
meeting the requirements given in the body of the 
paper and, in addition, having a pure wave-form output 
are necessary. Good amplitude and frequency control 
is especially important near resonances where the phase 
of the response is markedly affected by hunting. The 
Lockheed electric vibrators are satisfactory for ground 
testing but are too heavy and cumbersome to be car- 
ried in flight; air-jet-impulse or reaction-type exciters 
may provide the solution in that case. 

The prime need in instrumentation development is 
for an improved method of finding the phase of the 
forced response with respect to the applied force, since 
phase angle measurements made directly on sinusoidal 
traces of oscillograph records are both tedious and in- 
accurate. (A more or less successful attempt at solv- 
ing this problem has recently been made in Ger- 
many.’® 1”) Several possibilities suggest themselves 
immediately; there are doubtless many others. (1) 
A set of small cathode-ray oscilloscope tubes with 
simple single purpose circuits—one per pickup—from 
which the phase between the sinusoidal output of the 
pickup and that from a vibrator reference connection 
can be measured on an elliptical Lissajous figure. 
These ellipses, together with a frequency meter, might 
be recorded by a photo-observer camera, and the film 
could be read in a microfilm projector. (2) A pickup or 
galvanometer system that would record directly the 
real and imaginary components cf the response rela- 
tive to the exciting force. A wattmeter-type connec- 
tion between the vibrator and the measuring pickup™ 
may be used for this purpose. Advantages of this 
type of instrument would be its filtering characteristic 
to eliminate undesirable harmonics and the fact pointed 
out previously that the real and imaginary components 
of the total response are simply scalar summations of, 


* A similar situation may arise in the vicinity of the natural 
frequency of the elastic suspension of an engine; in a case pointed 
out by Myklestad,” nine out of a possible 17 wing resonances 
were calculated to occur in the region 0.9weng, < w < 1.5weng.. 
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respectively, the real and imaginary components of 
the responses in the various modes. (3) Means by 
which an electric signal having the frequency and phase 
of the applied force and a magnitude proportional to 
its amplitude could be added to, and/or subtracted 
from, the response of a pickup; thus, it would be pos- 
sible to obtain the amplitude and phase response by 
measuring the amplitudes of two sinusoidal traces, 
(4) Various instruments have recently been developed 
which indicate directly the phase angle between two 
electrical currents, and it should be possible to adapt 
these to oscillograph recording. The total amplitude 
of the component in the response wave form having 
the vibrator frequency would have to be simultaneously 
measurable. ‘ 


The test results give little reason to doubt the appli- 
cability of the theory presented, since the simpler cases 
check and all disagreements may be explained on the 
basis of the uncertainties already mentioned. Ob- 
viously, the assumption that is most open to question 
and may therefore require further study, is that the 
modes of vibration of an actual airplane structure have 
the properties enumerated at the beginning of this 
paper, which, as is well known, is rigorously true only 
for the normal modes of an undamped system. While 
the purely structural damping due to hysteresis in the 
material is believed to be small enough not to cause 
any major departures, this is probably not so for engine 
or control-surface modes, since considerable energy 
may be dissipated in the engine dampers or, respec- 
tively, in the friction in hinges or control cables. Theo- 
retically, it can be shown that the total forced vibration 
response of a linear mechanical system having large 
damping is also a vectorial sum of responses in various 
modes; moreover, the variation with frequency of each 
of these components can still be represented by the 
circles discussed in this paper and, consequently, iden- 
tical frequency response vector plots will be obtained. 
When large damping is present, the mutual interference 
between different modes increases, and, therefore, it 
becomes even more important to consider the phase as 
well as the amplitude response—i.e., to think in terms of 
vectors. A more thorough theoretical study of the ef- 
fect of damping will be made and several examples of 
experimentally obtained frequency response vector 
plots will be given in a later paper. 


The main purpose of ground vibration testing of air- 
planes is to check directly the natural frequencies and 
shapes of the various coupled modes calculated from 
mass and stiffness distributions. These modes are neces- 
sary when making flutter analyses or when it is desired 
to treat the airplane as an elastic body in the computa- 
tion of transient stresses such as those produced in 
landing or by air gusts in flight. The attempts made 
at verifying the accuracy of predicted flutter speeds or 
landing and gust loads through full-scale testing of air- 
planes have not been too successful in the past. In the 
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future, there can be little hope of obtaining more posi- 
tive results as long as it is not definitely proved that 
the actual modes of vibration of airplanes in still air 
can be satisfactorily calculated and, therefore, that 
proper account is taken of the internal structural forces. 
In the light of the work presented in this paper, the 
solution to the latter problem depends on the future de- 
velopment of more suitable vibration-measuring equip- 
ment with which to be able to identify and analyze 
airplane vibration modes accurately. 
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Letters to the Editor 


Dear Sir: 

Recent work for the American Iron and Steel Institute on 
light-gage steel sections has brought to my attention an article 
by George Gerard on “Effective Width of Elastically Supported 
Flat Plates” (JOURNAL OF THE AERONAUTICAL SCIENCES, Vol. 13, 
No. 10, October, 1946). 

On page 521 it is stated: ‘‘Since changes in half-wave length 
cause small variations in critical stress of plates with aspect ratios 
‘(a/b)’ greater than 2.5, it is probably permissible for practical 
considerations to establish 2.5 as the lower limit of applicability 
of the results presented.”’ 

In view of the work of A. J. Miles quoted by S. Timoshenko, 
Theory of Elastic Stability, the case of a plate rigidly supported 
along the two edges normal to the direction of compression, 
elastically supported along the other two edges (torsion of the 
supports neglected), it can be shown that there is an upper limit 
of aspect ratio beyond which the effective width falls off with extreme 
rapidity. Using the notation of Timoshenko, let Jp be the moment 
of inertia of one elastic support (including effective plate width 
in the determination thereof) and Jp the moment of inertia of the 
whole plate section of width 6. Then for our present purpose, 
the @ function is approximately (1 — y?)Ir/Ip. For long plates, 
m = 1 and the ¢ function @ = x(b/a) and is small. Then the 
v function is approximately y = Vé¢*. If y is set equal to 2z, 
then it can be shown that (a/b)o, the upper limit, is (a/b)o = 
1.25V (Ir/Ip)» somewhat approximately. For lesser values of 
(a/b), it can be shown that y varies nearly 6 per cent from con- 
stancy, as first one m value then another tends to govern. 


Hence, in this region, allowance must be made for such wave ef- 
fects in correcting data. 

Above (a/b), the effective width varies inversely as the square 
of the ratio (a/b)/(a/b)o, so that at twice (a/b)» the effective 
width is only 25 per cent of that for shorter lengths. This effect 
has not been noted heretofore, but has been observed on column 
tests. 

Reference to the JouRNAL, ‘‘An Investigation of Duralumin- 
ium Struts Under Compression,’’ by Warren S. Parr and Wallace 
M. Beakley, Vol. 3, September, 1935, shows that the theoretical 
plate failure line was departed from further and further with in- 
creasing b/w ratio of flange depth to channel depth (flange con- 
stant, hence depth decreasing and, for constant strut length, 
(a/b) increasing) starting at a fairly definite point. A check of 
the point of departure agreed excellently with the above theory 
of reduction in effective plate width with (a/b)» ratio in three of 
the four cases permitting a check, and the percentage decrease 
was of the right order of magnitude. It is of interest to note that, 
in the so-called region of torsional failure, neither the torsion line 
nor the above plate line (reduced for length) governed; rather it 
seemed that the failure was a mean of the two values. 

I have not made reference to E. Chwalla (Timoshenko, 
Theory of Elastic Stability, also Ingénieur-Archiv., Vol. 5, p. 54, 
1934) for the effect of torsion in the case of elastic support. Some 
slight variations have been noted. 

I have found that the tests of George Winter (‘‘Strength of 
Thin Steel Compression Flanges,”’ Proc. ASCE, Vol. 72, No. 2, 
February, 1946) on light columns consisting of two channels 
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without stiffened flanges, fastened back to back, exhibited the 
same effect. So far, the data available have permitted check of 
four cases in which the theory permits of accounting for the ob- 
served results. 

It should be remembered that, when m = 1 and (a/b) is greater 
than (a/b)o, the plate and elastic supports deflect as a single 
plate of the stiffness of the reduced section without supports along 
the sides parallel to the direction of compression. 

Investigation shows that many light columns in the so-called 
primary member category, so far as L/r ratio is concerned, may 
exhibit this effect and that some in the Euler range may be re- 
duced as much as 70 per cent in strength. Hence, the effect can 
prove to be important. 

The tests of Winter exhibit the effect of change from simple 
support at low values of (b/t) to clamped support at the higher 
values, much as suggested in Gerard’s Fig. 1. 

These remarks are purely preliminary to more adequate pub- 
lication but seem demanded by the occasional importance of the 
effect and the apparent failure heretofore to recognize it. 


Epwarp ADAMS RICHARDSON 
Publications 
Bethlehem Steel Company 





Dear Sir: 

I would like to add a few remarks to my previous letter con- 
cerning the So-Called Principle of Least Work Method. In con- 
sidering the application of the Theorem of Least Work to the 
solution of practical problems, two main cases must be distin- 
guished. The separation into these two cases was already being 
made by Péschl (Theodor Péschl, ‘Uber die Minimalprinzipe 
der Elastizitatstheorie,” Der Bauingenieur, Vol. 17, No. 17/18, 
p. 160 to 164, Mai, 1936. See also discussion of this paper by 
O. Domke, Der Bauingenieur, Vol. 17, No. 41/42, p. 459, 1936.) 
This distinction is also self-evident from the exposition of the 
definition of the So-Called Principle of Least Work Method given 
by me in my last letter, published in the October issue of the 
JOURNAL. 

The first case is when, in the system considered, we deal only 
with the concentrated, separate forces in the members of the 
structure. In this case the total strain energy of the system is 
given by a quadratic form which is always a positive definite one; 
therefore, the extremum sought must bea minimum. The author 
cited previously various attempts (found in several textbooks) to 
prove the existence of a minimum by pointing out that the sec- 
ond partial derivatives taken separately with respect to X, Y, Z, 
etc., are positive. This condition is necessary but not yet suf- 
ficient. The fact that the quadratic form is a positive definite 
is the only proof of existence of a minimum. (See: Theodor 
Péschl, Elementare Festigkeitslehre, p. 173; Berlin, Julius Springer, 
1936.) In this case the solution is given by equating to zero the 
first partial derivatives of the total strain energy with respect to 
the redundant forces. Péschl calls this case ‘‘Castigliano’s Prin- 
ciple.’’ This is a trivial case. 

The second case occurs when, in the system, there are continu- 
ous distributed stresses. The total strain energy is expressed in 
this case by single, double, or triple integral forms, and the 
equating to zero of the first partial derivatives of the total strain 
energy with respect to the unknown coefficients may not lead to 
a correct solution of the problem. Péschl recommends in this 
case as a constraining condition the equality of the strain energy 
and the work of the external forces (conservative system). 

Two methods may give the solution. The first is to assume the 
curve or the surface representing in a first approximation the dis- 
tribution of stresses as a two-dimensional curve function or a 
three-dimensional surface function of the system of coordinates. 
This choice may be made in many ways. The simplest may be 

















to base the distribution on experimental tests. The chosen 
functions possess certain unknown coefficients, which may be 
determined by the application of the Theorem of Minimum Strain 
Energy to obtain the proper solution. As an example, let us 
mention that this way was chosen by Prof. Eric Reissner in his 
solutions of the shear lag problem. (See: E. Reissner, ‘‘Least 
Work Solutions of Shear Lag Problems,’’ JOURNAL OF THE AERO- 
NAUTICAL SCIENCES, Vol. 8, No. 7, p. 284, May, 1941. This ap- 
proach has its full justification in the proposition put forward by 
the author in the previous letter. Before we apply the Minimum 
Strain Energy Principle we are assured that among the chosen 
set of functions there is one that gives a lower value of the energy 
than do the others. This method may be considered to be ap- 
proximate, but it is shorter. 

The other approach is to apply the methods of the calculus 
of variations. It is no longer a problem of seeking an ordinary 
extremum of a function. This well-known method is thoroughly 
explained in the literature (Courant, Hilbert, Sokolnikoff, 
Péschl, Joc. cit.). Strictly speaking, the use of the principles of 
calculus of variations to the solution of problems of this kind was 
recommended as far back as 1914 (Oskar Domke, “‘Uber Varia- 
tionsprinzipen in der Elastitatslehre nebst Anwendungen auf die 
technische Statik,’’ Zeitschrift fir Mathematik und Physik, p. 174, 
Mai, 1914). Thus, the general theoretical way of attacking a 
problem does exist. Strangely enough, although Courant-Hil- 
bert’s book was published in 1931, and Péschl’s paper in 1936, in 
many books on the theory of elasticity, theory of structures, 
strength of materials, aircraft structures, etc., published later 
than the above-mentioned works, there is nothing said about the 
vital and basic importance of distinguishing between these two 
cases. Also few books mention the necessity of the application 
of calculus of variations methods (see Sokolnikoff). There are 
many examples cited in which the continuous distributed stresses 
were found by equating to zero the first partial derivatives. But 
the enormous and insurmountable difficulties in finding exact solu- 
tions of differential equations occurring in variational calculus 
methods applied to such problems as shells, etc. (particularly with 
cutouts), were the reason for the author’s remark that the question 
of how to find such a function in any particular case is not yet 
solved in a rigorous way. Perhaps from an engineering stand- 
point it may be enough to apply some approximate solutions by 
use of fast computing machines (punch cards), although the com- 
putation difficulties seem to be enormous. 


+M. Z. KRZYWOBLOCKI 
* College of Engineering 
Associate Professor 
University of Illinois 





Dear Sir: ; 

I have two wind-tunnel items of interest; both concerned with 
blocking corrections. : 

(1) It was desired to test an NACA 0015 airfoil section 
through an angle of attack range of 180° in order to get forces and 
pressures for a helicopter blade, and some concern was felt about 
the blocking corrections when the 18-in. chord model was at an 
angle of attack of 90° in our-2'/2 by 9 ft. two-dimensional tunnel. 
Consequently, two additional models of 6- and 12-in. chord were 
constructed and tested. All models were mounted between flush 
wall discs. Leakage to the outside around the discs was prevented 
by water seals. The speed was set by a single piezometer ring 
in the tunnel settling chamber, referenced to the atmosphere. 
Although this method of speed setting is inferior to a two-ring 
system, it is believed that no serious effects resulted thereby, since 
pressures read on the 18-in. chord model and worked up, assuming 
clear jet g, agreed well with the force data. 

The force data as read included the tare and interference of the 


mounting discs. The true airfoil drag was obtained with a 
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Fic. 1. Drag data uncorrected for blocking of three 0015 airfoil models in 2'/: by 9 ft. two-dimensional tunnel. 
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Fic. 2. Drag coefficient [= drag/(d*/4)q] of flat discs in open test section, uncorrected for blocking. 
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wake survey rake at zero lift, and the tare and interference thus 
obtained were subtracted from allreadings. This admittedly neg- 
lects the portion of the tare and interference that varies with 
angle of attack. The error, particularly in the 90° angle-of- 
attack range, is believed to be a small fraction of the total drag. 
The results, uncorrected for any blocking whatsoever, are shown 
in Fig. 1. Although reference 1 indicates that a wake blocking of 
44 per cent qg, and reference 2 indicates that a wake blocking of 
13.0 per cent g, could arise, apparently no appreciable nlocking 
exists, or the three widely different models would surely show it. 

(2) When testing rotors in the vertical descent condition in 
an open throat wind tunnel, a flow pattern not greatly dissimilar 
from that about a round solid disc is developed. Reference 3 
found a blocking effect of 34.5 per cent g for a disc of 0.5D. In 
order to check that figure, discs from 0.11D to 0.555D were 
tested in the Georgia Tech Nine Foot open throat wind tunnel. 

The discs were mounted at the tunnel centerline 5 ft. from the 
entrance cone with a single shielded support that came up well 
behind the disc. The tare of the support (which proved negligibly 
small) was evaluated by supporting the discs externally with a 
small clearance about the center mount. The interference of the 
support behind the disc was believed to be negligible from experi- 
ence gained testing rotors with supports in their wake. 

The tunnel air speed was set using a single piezometer ring in 
the settling chamber, referenced to the atmosphere. 

The results are shown in Fig. 2, uncorrected for any wall or 
blocking effect. Although a considerable scatter exists for the 
10 and 20 m.p.h. values, doubtless due to the difficulties of hold- 
ing the tunnel speed steady at such extremely small values of the 
dynamic pressure, the evidence is that the disc drag is unaffected 
by the tunnel boundaries at least up to discs of 0.555D. 

It is in order to explain the discrepancy between these tests and 
the theory. Briefly, the preponderance of blocking for a flat plate 
is the ‘‘wake blocking’’ type. The theory is based on the true 
assumption that the wake behind a model is at a lower speed 
than free air} Consequently, from continuity considerations, the 
air surrounding the wake is at a higher velocity and lowered 
pressure. It is then supposed that the lowered pressure reacts 
back through the wake to the rear of the model and tends to draw 
it downstream, increasing the drag. Glauert* pointed out that 
the entire stream is free to contract in an open jet, and, hence, 
no wake blocking should appear. This was.sibstantiated in 
Item (2). However, the author feels that far more mixing must 
take place in the closed jet than was originally considered, and, 
hence, the lowered pressure that was assumed is not actually 
realized. Hence, wake blocking is, according to these data, an 
extremely small factor. 

It is planned to present more complete data on the tests of 
Item (1) at a later date. 


REFERENCES 


1 Glauert, H., Wind Tunnel Interferenice on Wings, Bodies, and Airscrews, 
R & M 1566, 1933. 

2 Thom, A., Blockage Corrections in a‘Closed High Speed Tunnel, R & M 
2033, 1943, 

* Wheatly, John B., and Bioletti, Carlton, Wind Tunnel Tests of a 10 Foot 
Diameter Gyroplane Rotor, TR 536, page 7, 1935. 


ALAN PopE 
Associate Professor 
Daniel Guggenheim ScHool of Aeronautics 
« Georgia School of Technology 





Dear Sir: 

Some of your readers will probably recall séeing an article in 4 
popular magazine concerning the abilities of a bumblebee to sus- 
tain himself in flight and presenting in particular some confusing 
formulas purporting to prove the author’s contentions. The 
article, which was later reprinted in several places, was supposed 








excessive wing loading. The old anecdote quotes an ‘‘aerodynami- 
cist”’ to the effect that W/S is too great for flight but goes on to 
say that the bee, ignorant of. aerodynamics, flies nevertheless. 

The magazine story, as I have mentioned, was supposed to kil] 
this old jibe at aerodynamicists by presenting experimentally 
deduced formulas that would show how the little insect might 
produce adequate lift in spite of this supposedly large unit load- 
ing. Your readers may be interested to know that the formulas 
actually had their origin in a serious experimental investigation. 
Unfortunately, they were somewhat distorted by a typesetter 
and, when they appeared in print, were almost certainly dismissed 
as nonsense by any engineer who read the story. 

The experimental investigation was made by E. von Holst 
at Géttingen, as reported in Luftfahriforschung of 1943.* In this 
remarkable series of tests, he placed in an air stream an elliptic 
wing of 3.2 cm. span and caused it to perform simultaneous beat- 
ing and pitching oscillations. Measurements were made of lift, 
drag (or thrust), and power required, all as functions of the mean 
angle of attack a, the angle y between the windstream and the 
axis of oscillations, the pitching amplitude x, and the reduced 
frequency or advance ratio \ = v/u (where u = mean tip speed 
and v = stream speed). The beat amplitude ¢ was constant. 

To assist the author of the magazine article, I selected a typi- 
cal case: 

a = 45°, y = 45°, ¢ = 68° 
for which, the experimental data show, the lift coefficient is ap- 
proximately represented by the straight line 


0.06 + 0.17X & ka [= lift/(p/2)u?F] 


x = 40°, 


(p = air density; F = area swept in one beat = 2ns*(y/2r); s = 
span). When this expression is used to calcilate the lift, one 
finds 

lift = (p/2) F|0.06(2sng)? + 0.17 (2sn¢)v] 


where n is the frequency of the oscillation. 

It is clear that the bumblebee can produce adequate lift at 
zero flying speed by simply flapping his wings fast enough (as 
must appear obvious). Moreover, he can reduce as he picks 
up speed. von Holst’s results include data from which the thrust 
and power required can be evaluated. For the particular case 
mentioned above, for example, positive thrust is obtained over a 
considerable range of the parameters. 

With such complete data at hand, there appears to be no need 
for basing any rash conclusions on the value of any one aero- 
dynamic parameter! May we consider the ancient libel on aero- 
dynamicists refuted forever? 
~ * von Holst, E., Messungen an kleinen schwingenden Flachen (Gréssenord- 
nung des Insektenflugelns), Luftfahrtforschung 20, Lfg. 5, pp. 137-146, 1943. 

W. R. SEARS 
College of Engineering 
Cornell University 





Dear Sir: 

I am prompted to send a letter to your JOURNAL in.answer to 
the letters of H. Reese Ivey (page 450) and Prof. Ralph H. Upson 
(page 469) in the August, 1947, issue of the JOURNAL. 

Whatever quibble I might have regarding the value of propul- 
sive efficiency , asa performance criterion is not based on the 
fact that np fails to be a plausible and rational quantity. ‘It does, 
in fact, serve to determine the ratio between missile energy and 
energy from fuel expenditure, as it is obviously meant to do. I 
am not inclined to quibble about the plausibility of the deriva- 
tion, as between the forms of n, proposed by Professor Szczeniowo- 
ski and those of Mr. Ivey and Professor Upson. 

My only comment is that neither form of n, serves to answer 
any pertinent question in Aeronautical Engineering. This is to 
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LETTERS TO 


say that economy of energy expenditure per se is not a practical 
consideration in the design of missiles. The primary question 
that is asked in our scieace is: How fast can a given pay load be 
carried how far? In answering that question the aim is to reduce 
the mass ratio m,/ m to as small a value as possible or, for given 
m,/'m, to sead the net mass as fast and as far as possible. Whether 
in so doing the energy expended is much or little is of no more real 
jnterest than how much the fuel costs. Such mundane consider- 
ations as are used in deriving », are of primary importance in 
stationary power-plant engineering but not in Aeronautical 
Engineering. 

I recognize that in the preceding statements I am simply voic- 
ing an opinion, and there may be some who disagree. There is 
no reason why the course of the energy transformationssshould 
not be traced, if anyone is interested, whether such an activity 
helps design rockets or not. On the other hand, no one should 
deceive himself: if there is a criterion like n, that can be large or 
small or negative, that is decreased by an obviously advantageous 
change in design of rocket, that can be increased by the same 
change at a different missile speed, then such a performance 
criterion is suspect if not useless. The fact that a similarly devel- 
oped criterion is useful in other engineering fields is not a good 
enough reason for retaining it. 

Davin T. WILLIAMS 
Associate Professor 
Department of Aeronautical Engineering 
Univeisity of Michigan 





Dear Sir: 

The publication of Mr. Wang’s article, ‘‘Elastic and Plastic 
Bending of Beams’’ (JOURNAL OF THE AERONAUTICAL SCIENCES, 
Vol. 14, No. 7, p. 422, July, 1947), raises a question of editorial 
policy which I believe is of interest to all scientists and engineers. 
This concerns the introduction of new words in scientific or 
engineering terminology. A few simple rules are suggested, as 
follows: 


(1) The word should not be based on the name of the com- 
pany in which the author happens to be employed. 

(2) The need for a new word should be clearly shown—i.e., 
the possibility of using existing words, or combinations of them, 
should be thoroughly investigated. 

(3) If introduced at all, the new word, or words, should be at 
least partially self-explanatory. 

(4) The physical significance of the new term should be 
clearly explained. 


Applying these suggested rules to Mr. Wang’s article we find 
the following: 


(1) The word “‘Norstress”’ is derived from Northrop Aircraft, 
Inc., where the author was working at the time the term was 
originated. 

(2) The new term was introduced “for convenience.” It 
may represent a convenience to those who originated it and use it 
in routine work, but it is an inconvenience to everyone else who 
reads the paper. 

(3) Such terms as “equivalent bending stress” or “reduced 
bending stress’ would be preferable, for reasons shown below. 

(4) The only definition given by the author is a ‘‘term which 
represents a measurement of maximum normal stress (remote 
fiber stress) of cross sections of beams due to bending.’’ The 
actual physical significance of the term in one of its forms 
(M/BH?) can be revealed by writing this expression as (M/H)/- 
(BH). The term M/H represents in equivalent axial load, P, 
which is simply the bending moment divided by the maximum 
depth of the beam. The term BH represents the cross-sectional 
area, A, of a rectangular beam having the same depth and maxi- 
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mum width as the section in question. The so-called ‘‘Norstress,”’ 
in this case, is actually the fictitious axial stress obtained by con- 
verting the bending moment into a couple composed of axial 
forces that are assumed to act at the top and bottom of the cross 
section. The axial force is then assumed to act uniformly over 
the entire area of the superimposed rectangle. If the actual 
cross section is a rectangle, the ‘‘Norstress”’ turns out to be one- 
sixth of the maximum fiber stress obtained from the conventional 
bending theory for the elastic range. This suggests that the term 
(M/BH?) might well have been multiplied by 6, in which case the 
quantity could be defined as ‘‘the equivalent elastic bending stress 
in the superimposed rectangular cross section.” 

Lack of a clear physical basis for the new term has unfortu- 
nately permitted the author to introduce the term (./4/R*) for 
circular and tubular beams. To be consistent, the diameter, D, 
should have been used; in fact, the term (\4/BH?) would have 
served just as well. 

As a result of this, the ‘‘ Norstress”’ for sections of circular shape 
is eight times that for any other cross section having the same 
overall dimensions! (I wonder how the author would have 
treated an elliptical section.) The mere fact that this situation 
was not corrected before publication should be proof enough that 
the introduction of arbitrary concepts of this type is often likely 
to produce confusion rather than clarity. 


The whole thing could be cleared up by establishing a ratio (1 
won’t attempt to coin a word for it) between the actual stress at 
the outer fiber and the corresponding stress in the outer fiber of a 
member having a rectangular cross section which just encloses the 
actual cross section. This ratio would be similar to the present 
concept of a ‘‘form factor,” except that it would include both the 
inelastic effects and the geometrical characteristics of the cross 
section. As such, it would be convenient in stress calculations, 
because it would eliminate the procedure of calculating the elastic 
stress and then correcting this to account for inelastic effects. 


F. R. SHANLEY 
Consulting Structural Engineer 





Dear Sir: 

The numerical procedure described by Marvin Stern in his 
article in the August, 1947, issue of the JOURNAL OF THE AERONAU- 
TICAL SCIENCES (page 481) is a valuable contribution toward the 
solution of the ‘‘one-directional”’ shear lag problem and, with the 
slight extension mentioned at the end of the article, the shear lag 
problem involving lateral movement. 

The procedure suffers from some lack of generality, however, 
inasmuch as it cannot be applied directly to problems in which 
some of the boundary conditions are given in terms of displace- 
ments rather than stresses. In many cases this lack of generality 
does not represent any practical limitation on the applicability 
of the method, since a boundary having zero displacement can 
often be regarded as an axis of symmetry of a hypothetical sym- 
metrical structure and thus eliminated from consideration as a 
boundary. In some cases, however, particularly in problems 
involving skew boundaries or problems in which the prescribed 
boundary displacements are not zero, such an artifice is not pos- 
sible and some other way has to be found of handling the displace- 
ment boundary conditions. 

I should therefore like to describe a slight modification or exten- 
sion of Mr. Stern’s approach that makes possible the explicit 
handling of displacement, as well as stress boundary conditions. 
This modification involves the use of displacement equations 
rather than stress equations and the conversion of stress boundary 
conditions into conditions on rate of change of displacement. 

The equation of equilibrium for an element dx of stringer n at 
lattice point (m, m) is first written in terms of u-displacements 
as 
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o*u Un-1, m — Unm Uns mw — Un+ism\ 
Aut( >) Ti ue r r )=0 
(1) 


If the u-displacements in the mth stringer are assumed to vary 





according to a second-degree parabola between lattice points . 


(nm, m—1) and (n, m+1), then 
(07u/Ox*) ny m 
can be expressed in terms Of tn, m—1, Un, m, ANd Un, m41 2S 
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(2) = Zar mee +1(#) | eo 


(The assumption of a parabolic variation is also implied in the 
equation above Fig. 1 in Mr. Stern’s article.) Substitution of 
Eqs. (2a), (2b), and (2c), respectively, into Eq. (1) and solution 
of the latter equation in each case for un, m give the following three 
equations for tn, m: 


Un, m2(1+ 8?) = Un, m-1 + Un, m+1 + Bun-1, m + Buns, m 


re) 3 
2) ni + Un, mai + 2 Bun-1, m + 


(3a) 
Un, m(1 +36?) = -1( 


(3b) 


3 ty 
2 n+l) m 
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Ou 3 
Un, m (1+38?) = Un, m—-1 +1 (=) + - BPun—1, at 
Ox ny m+1 2 


: Bt n+1, m (3e) 
(If n is an edge stringer, Eq. (1) and therefore Eqs. (3a)—(3c) must 
be changed somewhat to eliminate the »—1, m or n+1, m 
term.) 

Eq. (8a), which relates the u-displacement at lattice point 
(n, m) to the u-displacements at four neighboring lattice points, 
takes the place of Eq. (3) in Mr. Stern’s article. If either (, m—1) 
(n,m + 1) isastringer-end point at which the stress os rather than 
the u-displacement is known, the known stress can be converted 
into a Known value of 0u/dx by use of Hooke’s law, and Eq. (3b) or 
(8c), respectively, must be used instead of Eq. (3a) at lattice point 
(n, m). 

By writing Eq. (8a), (3b), or (3c) at every lattice point, a 
set of equations is obtained that involves as many unknown 4- 
displacements as there are equations. These equations can be 
solved for the u-displacements by means of the iteration procedure 
described by Mr. Stern. It is only necessary to assume values for 
the un, m’s and to adjust each value using Eq. (3a), (3b), or (3c), 
depending on what conditions are known at the neighboring lat- 
tice points. After the displacements have been determined, the 
stringer stresses can be calculated from the rate of change of dis- 
placement, which is the strain. 


CHARLES LIBOVE 
Structural Engineer 
National Advisory Committee for Aeronautics 








Sciences, Vol. 14, No. 6, pp. 351-358, June, 1947. 
On page 352, Eq. (12) should read as follows: 





Errata 


The following corrections to his article were submitted by John W. Miles, the author of the paper ‘‘The Aero- 
dynamic Forces on an Oscillating Airfoil at Supersonic Speeds,” which appeared in JOURNAL OF THE AERONAUTICAL 


G(x, &) = eM O-8) Ty [x(x — &)) 
On page 357 in Fig. 9, the ordinate should be: —k*J,a. 


The real and imaginary parts of the flutter derivatives plotted in the figures are referred to the quarter chord, 
for M = Oand M = 0.7; accordingly, direct comparison with the supersonic results in the figures is possible only 
in. the case of R,, and J,,, the rotation @ and the stalling moment being about the quarter chord. 
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The Linearized Subsonic and Supersonic Flow 


About Inclined Slender Bodies of Revolution 


E. V. LAITONE* 


Cornell Aeronautical Laboratory 


SUMMARY 


Rigorous first order linearized solutions are obtained for the 
subsonic and supersonic flow about bodies of revolution. The 
surface pressures are found by means of a Taylor’s Series ex- 
pansion in terms of the cross-sectional area. The resulting simple 
expressions are shown to be in satisfactory agreement with exist- 
ing experimental data. It is proved that previous solutions of 
the linearized supersonic flow equations have retained terms 
that at high Mach Numbers are incompatible with the basic 
potential equation of compressible flow. 

Extremely simple relations are developed for the lift and mo- 
ment of bodies of revolution in inclined and curvilinear subsonic 
and supersonic flow. These relations are in good agreement with 
experimental measurements. The commonly used linearized 
solution for the supersonic flow about an inclined body of revolu- 
tion is incorrect, since nearly all of the terms retained are of a 
higher order than the terms omitted in neglecting the effect of the 
velocity component due to the axial flow and in not considering 
the complete boundary conditions. 

The rigorous analysis indicates that the linearized theory can- 
not correctly predict the variation of the surface pressure coef- 
ficient at high Mach Numbers. 


SYMBOLS 


a = local velocity of sound (ft. per sec.) 

Mach Number = Uo? + V0o2/a A ta/a 

VM? — 1 

V1 — M 

isentropic gas constant = 7/5 for dry air 

= local velocity along x axis (ft. per sec.) 

= local velocity along r axis (for w = 0) (ft. 

sec.) 
Uo, Vo = free-stream velocity conditions at an infinite dis- 
tance upstream of the body (ft. per sec.) 

free-stream angle of attack, positive nose up = 
tan“ (ve /uUo) AR Vvo/te (rad.) 

= pitching velocity, positive nose up (rad. per sec.) 

local pressure coefficient = [(p — po)/q] 

= local static pressure (Ibs. per sq. ft.) 

free-stream static pressure (lbs. per sq.ft.) 

= p(ta* + V0*)/2 A pliw?/2 

= mass density (lb.sec.? per ft.*) 

(x, r,w) = cylindrical coordinates for body of revolution, 
Fig. 1 

cross-sectional area of body of revolution = 
mr?(x) (sq.ft.) 

= dS(x)/dx = 2xrr’ (ft.) 

dr(x) /dx 

length of body of revolution (ft.) 

= base area of body of revolution, S for x = L 

(sq.ft.) 
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* Principal Aeronautical Engineer. 
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cone or nose semivertex angle (rad.) 

resultant transverse force per unit axial length 
(Ibs. per ft.) 

= lift coefficient = lift/gA 

pitching moment coefficient about nose (x = 0) = 
moment/gAL 


6 = 


oo 
i 


PO 
ui 


INTRODUCTION 


ire PRESENT DEVELOPMENT was started as an at- 
tempt to justify the linearized supersonic flow solu- 
tions for a body of revolution as introduced by Karman- 
Moore! and Tsien.2 A mathematically rigorous de- 
velopment of the linearized supersonic flow solution was 
attempted in order to indicate its limitations and to ex- 
plain some of the anomalous behavior of the Karman- 
Moore and Tsien solutions. 

The analysis shows that in general the effect of high 
Mach Numbers cannot be predicted by the linearized 
theory for the supersonic flow about a sharp-nose 
slender body of revolution. The difficulties encoun- 
tered through the use of the Karman-Moore or Tsien 
linearized solutions are shown to be due to the fact that 
they implicitly contain terms that are essentially of an 
order higher than those that have been already neg- 
lected in the process of linearizing the basic equation 
of a compressible potential flow. These terms contain 
the Mach Number in such a manner that the numerical 
value of these terms becomes larger than that of the 
lower order terms as the Mach Number increases. 
The mathematical development indicates that in gen- 
eral the additional work required by the Karman- 
Moore and Tsien methods is not justified. In fact, 
the additional work involved would generally lead to 
erroneous results, at least for Mach Numbers greater 
than +/2. 

The analysis presented here was carried out during 
1946 at the Cornell Aeronautical Laboratory as a part 
of the investigation of the dynamic stability of missiles 
at subsonic and supersonic speeds. The work was 
done under the sponsorship of the Air Matériel Com- 
mand, Wright Field. 


THE SUPERSONIC AXIAL FLow AsBouT A ‘Bopy OF 
REVOLUTION 


The potential equation for the axially symmetric 
flow of a compressible fluid in cylindrical coordinates 
(Fig. 1) is* 
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Tweee DIMENSIONAL FLOW where ¢ is a perturbation potential superimposed on a 
w/ r uniform constant flow field, Eq. (1) becomes 
2 
Un >a eo. x ei ¥ [¢=(1 sal sues 
«x lan’ 2 
Uo z Tr 
Un pan “A> E (2 + $=) + PR . 
Op a a - (3) 
Uo 2 \Pr 
s(t #8] 
(er,r) (x1) a ae 
a= 6.* ings ly _ 1)/2}(2u.¢z + o:” + ¢,”) 
mi lan™ i t 2 
(20) C a DD re which is linearized to the first order as 
— uy ee 2 
On THE MACH CONE (Se) $ Prr + (1/r)¢, + (1 M baz 0 + 0(¢ ) (4) 
Fic. 1. where for supersonic flow we define 
&,? &, ) = (u."/a.7) -1= M*-—-1>0 
6, (1- %) + 4,1 - +2 - 
a r K4rman-Moore! used the following as a solution of 
®,, "Eq. (4): 
26,, ro =Q} (1) q ( ) 
@ pe ee See cosh ~!x/ Br 
¢ = —(1/4r) f(x — Br cosh z)dz (5) 
y¥-1 1 
a* = a, 5) (u.*— ,? — $,? — 5%) or with 
Then if ® is taken as § = x — Brcoshz 6) 
6 
b=u.x+¢ (2) $ = —(1/4x) Lt” f(e)dt/V (x — £)? — BY? 





Eq. (6) is an improper integral, since the denominator vanishes for the upper limit, and, consequently, all pre- 
vious differentiation of the perturbation potential has been accomplished by operating on the form given in Eq. 
(5). However, some new useful relations can be obtained by differentiating Eq. (6) in the following manner: 


















































‘ &  f()(E — x)dé f(ii)(dii/dx) 
—4r¢, = L as, 
Y vim, | 0 [(x oF £)? vex B*r?}8 f V (x — £,)? — Bs 
¥ * [f(é) — f(&)) f(&) f(&) f(&) (dii/dx) \ 
L d ‘a 
um 1, Ge — e)? — Bra F008 + ee Va b)- BYt | VG — by BY 
=—Br (f(t) — f(x — Br)] S(« — Br) 
igo a Bepa tt 
ee ee fede flé)(d/dr) 2 
ht, taal 12 rf [je — 6? — Bri + 7 gy? — Bye 
ma ay, (2 © =sde_—__e/ fl) Mi = AKG) SB) dV ar “ 
&i—>x — Br [(x — &)? — B%}8/2 ~ y/y2 — Bey? ~ Ve — &)? — Bet | V(x — &)? — BY? 





(8) 





py [7% UO fle - Bde (x/nfle — Br) 
"Jo =e — 8? — Bet” V2 — Be 
That this is a valid mathematical procedure is easily shown by direct substitution into Eq. (4). Of course, Eqs. 
(7) and (8) could be obtained by differentiating Eq. (5) to get ¢, and ¢, and then integrating these by parts to get 
Eqs. (7) and (8), respectively. 
Now if f(x) is a continuous function, 


a 
A) — fle — Br) = YEO) OPM poy 


n=1 





and Eqs. (7) and (8) become 
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—4r¢; = {[ 108 (—_ =) = we = Fe) 9 (x) — [Vx — BY Br), x) 


- 3 BY? Br B*r?(x — Br) x a ee 2 i ey 
hs wea) eee pinching 
shila [2 ( Br ) 2x5 + 5x*B*r? — 15xB*tr4 + er 
9 R2y2 2 _ R22 4y4 ————— x 
2B Vx cil Fe 24 + 8 a x — Wx? — B%? 8V x? — BY? 
(5) - B B? 2\ —1/2 


120 x? 




















- {5[ |r oui E Br yee 
vith = ELS ee i Irs Bll) + Br) oe Veni 2 | 


Br(x — Br) |f’’’(x) 3 Br gree oy 
‘< id Bl 3 Br tog (= Taam) + 5Va — Bet oO - 


. / 
oid I ihc Vx? — By? 

- — BY (x — 2 Fw _ f(x — Br) BYy?\-1/2 
B*r f (x? — B*r?)*/2 + 2B*r? Vx? — BY? — Vee ae ane 720 Bh Trig ero 1— >) (10) 








Then, expanding in powers of r, 


ii 1 1 , )— vt xP _ ) one x$ (O( ) x (5) 


724... + 0@%| an 
or = rahe wy +(F log =) ot ae ober | (12) 
The boundary condition for the meridian line is ¢, = ur’ + ¢,r’; therefore, to the first order, 
f(x) = 4xu, rr’ = 2u.S’ (13) 
where 
S = xr? (14) 


Then, for the second order, 


2y2 


B 
fle) = dar" un + 45) — (FF tog =) 77%) + OB) - 


2 
, Br\ (3B?rr"’ — B*r*r’”’ - v. 
aS log °F oo vw Or! [i tr ) + O(r?) (15) 





Therefore, to the order r‘ log r = 0, Eqs. (5) and (6) become 


Au 1 cosh ~1x/Br ws __ S(O) _ =|- qf SI! (b)dé j 
ite i S''(x — Br cosh z)dz + Vat — B on Va? shy Be + 


__ Ss’) 1] - pip: ~ Pr [S'(é) — Si — Br)] S@ — Br)\ ., 
Ve Bel ~ lJ, [e— 6? — Beppe F— E+ 7 St (16) 














Av 1 cosh ~1x/ Br s” B ‘. _— (x/ nS: (0) 
Tig 518 f (x — Br cosh z) cos z 7 oe 2 
ar! =—Br S'"(H(E — x)dt_ ——_(x/r)S so o 
2a 0 V(x - — ¢)? B?r? Vx! 2 — BY Ass. 
1 *~Pr [S'(&) — S'(x — Br)l\d& — (x/r)S'(x — Br) 
>) ee ee EY ees Po ere - 
ros f [(x — &)? — B%y2}3/2 V x2 — By? ; (17) 


Similarly Eqs. (11) and (12) become 


4 


Au 1 ae Br ”? x sr x? Ss x? (5) ad So 3 > 
hie oe 1 + dog 5 JS” = gS" + 558° — gS? + eS + ---- +007), 18) 


pe 2Qr 
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a= -3(2)-()+( 


which, in view of Eqs. (18) and (19), becomes 
C, = —2 (Au/u,) 
or 


i 


5’ Br x x? x x4 
PS a ge WS eens © oe. +, ts = ie 
G 142 (1 + log : og 2)s 9 Si + 12° SO + Ss? + .. 


Fig. 2 shows that Eq. (20) gives a solution almost 
identical to the one obtained by K4rmaén-Moore.' 
However, the terms of order r? contain the Mach Num- 
ber and cannot be retained at high Mach Numbers 
unless there is an additional justification for using 
terms of an order neglected in obtaining Eq. (4). A 
comparison with the experimental data for an infinite 
cone shows that at high Mach Numbers only the first 
order term may be retained. Eggs. (5) and (15) yield 
¢ 1 


Uno 2a Jo 


cosh ~!x/ Br 
S’(x — Br cosh z) dz + 


O(rtlogr) (21) 


Then, for an infinite cone with a semiwedge angle of @, 
6 = r/x, S = 26x? 


Eq. (21) becomes 


ka = —6 (« cosh! = — Vx? Br) + 0(64 log @) 
r 


Ue 
(22) 


Substituting Eq. (22) into Eq. (4) shows that it is 
satisfied if terms of the order 6‘ are neglected. This is 
in agreement with the approximation for f(x) in Eq. 
(15) and with the work of Lighthill.4 However, Eq. 
(3) is satisfied by Eq. (22) only if terms of the order 
6? are neglected; consequently, the linearized solution 
for the infinite cone is valid only to the order 6? = 0. 
Therefore, for an infinite cone, Eq. (20) must be re- 
duced from 


C, = —20? log (B@/2) — 0? + 0(6* log 8) 
to 
C, = —26? log 6 + 0(6?) 


The neglected term — 26? log (B/2) is of order 6? but is 
also a function of the Mach Number in such a manner 
that as M increases this neglected second order term 
can become numerically larger than the first order 
term for a finite 0. Therefore, without further justi- 
fication, only the first order term S”/m log (r/x) is a 


A 
+. = r'(x) + O(r? log r) 


The pressure coefficient for compressible flow is given by 





1947 





=) (M? — 1) + 0(Au', Av‘) 


— (r’)? + O(r4 log r) 


72 480 


|- (r’)? + (73) ' (20) 
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r — ee 
org i im = — 
. a SA 
o| fa 20a Pes" = 
-.] f 4 
| Pa 
-.2 — EXPERIMENTAL DATA 
| (REF 1) —~S 
» 1 1 
° 1 2 3 4 5 6 7 
X INCHES 
Fic. 2. 
valid representation of the surface pressures. By di- 


rect substitution into Eq. (3), the exact potential equa- 
tion for compressible flow, it may be shown that 
S”"/m log (Br/x) is also a valid solution as the Mach 
Number approaches unity. It is found, however, that 
the additional term log B cannot be retained in a first 
order approximation as the Mach Number approaches 
and exceeds +/2, not even if the Mach number ap- 
proaches infinity. The rigorous linearized first order 
solution for the surface pressures on an extremely 
slender sharp nase body of revolution in supersonic 
flow is therefore given by 


oa Br 


7 -—-—ie—, l1<M<vV2 
ra x 
sr (20a) 
r 
C, =~ — — log-, M> 2 
va x 


In Fig. 2, Eq. (20a) is compared with Eq. (20) and the 
experimental data from reference 1. Fig. 3 compares 
Eqs. (20) and (20a) with the Karman-Moore! solution 
and the more exact Taylor-Maccoll*: solution for an 
infinite cone. It is seen that, at least for the cone, Eq. 
(20a) is the best of the linearized solutions. For the 
infinite cone the Karm4n-Moore! solution shown in 
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Fie. 3. 


Fig. 3 is that given by reference 6 as 
C 26? cosh~! (1/Bé) 
2 V1 — B*™? + 6? cosh “1(1/B8) 


It should be noticed that this expression fails to give a 
solution when B — 1/6 as indicated in Fig. 3. Also 
shown in Fig. 3 is the limiting asymptotic value of the 
pressure coefficient due to the conical adiabatic shock 
wave of a very slender conical nose when the Mach 
Number increases without limit as given by the author 
in‘ reference 7: 


C, > (vy + 1)? = (12/5)6 (23) 
It is seen that the rate of decrease of pressure coef- 
ficient as given by the K4rman-Moore' solution is un- 


doubtedly excessive. 
The wave drag is given very simply by Eq. (20a) as 


drag force 1 L , ) 
= — = APS’ (x)dx 
* gA i Me 





1 L Br 
ne ae , ” —dx, 
eek S'(x)S (x) log — 0 | (24) 
1< M< v2 


L 
=~ — 9 I S’(x)S" (x) log “dx, M>~/2 
TA J, x } 





This gives less wave drag than the relation given by 
von K4rman,® 


1 ad £ Ss" A ad 1 
o-saf f (x)S’"(8) log 


However, if Eq. (20) is substituted into 
Cp = (1/A) fo"C,S’ (x)dx 


the wave drag coefficient will be found to be in agree- 
ment with that given by Eq. (25) for bodies of revolu- 
tion whose cross-sectional area S is a uniformly con- 
tinuous function, at least for Mach Numbers near +/2 
or at any Mach Number for sharp nose bodies that 
have either r(L) or r’(L) equal zero. This is so because 








Lede dxdt (25) 
x—& 


(25a) 


Eqs. (20) and (25a) show that the wave drag coefficient 
given by Eq. (20) is also independent of Mach Number 
if S’(L) is equal to S’(O). For a sharp nose body 
S’(O) is zero; therefore, according to Eq. (20), S’(L) 
must also be zero for the wave drag coefficient to be in- 
dependent of Mach Number. 

Although the rigorous mathematical analysis can 
justify only the first order solution given by Eq. (20a), 
it still may be possible to justify the use of Eq. (20) at 
least for values of B not too much greater than unity, 
since the terms of order r? only contain the Mach Num- 
ber effect in the form log B, which does not increase 
extremely rapidly with Mach Number. It is easily 
shown by comparison that, for any body of revolution 
whose cross-sectional area is defined by a continuous 
function of x, Eq. (20) will give the same surface pres- 
sures as given by the von Karman” * method for B = 1. 
The von K4rman solution, however, will overestimate 
the effect of an increasing Mach Number, since the 
terms of order r* which are retained in it are all of the 
form B’rS” [see Eq. (11)]. The addition of the terms 
of order M?r* greatly overestimate the effect of Mach 
Number as B exceeds unity. This may easily be shown 
by comparison with experimental data or the solution 
obtained by the characteristic method of Sauer.* 


Lighthill* showed that Eqs. (16) and (19) could give 
a solution of Eq. (3) to the order r‘ log r = 0 in the 
following form: - 


1 x — Br 


ie S'(ode 
> Jo V(x — §)? — B%? 





— [r’(x)]? (20b) 





However, the present development proves that Eq. 
(20b) still contains terms of the form M’*r* which 
greatly overestimate the effect of Mach Number so that 
Lighthill’s‘ results are also applicable only for M <+/2. 

If in the supersonic flow linearized solution terms of 
the order r? may be retained, at least for bodies other 
than an infinite cone, then Eq. (20) may be used for 
the surface pressures if the free-stream Mach Number 
is greater than M,,,,. for a detached shock wave and if 
S(x) is continuous. If there are discontinuities in the 
derivatives of the cross-sectional area, the von Karman! 
solution may be used with B = 1 (M = +/2). Eqs. 
(16) and (17) may be used to find the velocities any- 
where in the flow field behind the Mach cone produced 
by the sharp leading edge of a body of revolution. For 
example, it will be shown later than an integration of 
Eq. (16) indicates that directly behind the body there 
exists a wake with a positive pressure increment similar 
to that existing in subsonic flow. In a similar manner, 
Eq. (16) may be used to determine the wind-tunnel wall 
and other interference corrections where the values of 
r are extremely large. 

It is important to note that all of these supersonic flow 
solutions are applicable only if the flow is everywhere 
supersonic. These solutions will not apply if there is a 
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detached nose shock wave. However, the minimum quently, the detached nose shock wave will seldom 
Mach Number for a detached shock wave in front of a occur in the final steady state supersonic flight of a 
conical nose with a vertex angle of 20° is 1.04;5 conse- properly designed sharp nose body of revolution. 


THE SUPERSONIC FLOW ABOUT AN INCLINED Bopy OF REVOLUTION 


Without axial symmetry the potential equation for the motion of a compressible fluid in cylindrical coordinates 
(Fig. 1) is* 


( ©,” ,? a ®,,” ®,®,, ®,®,, 
0 1 Gal It She — 26,,—— — 26, — — 
a r°a? ra? ra? 


am r 











©,h, , ®,,? —1 1 
Wh, ++ “(1 + =) =0, a? =a,?+ eC +-v.? — %,? — 6, — - #2) (26) 
a 2 7 
Then, for the linearized mathematical procedure, ® is taken: 
® = (ux + v.27 cos w) + O(x, 7, w) = Ua (x + racosw) + ¢ + O(ra’*) 
where ¢ is a small perturbation potential superimposed on a constant uniform slightly inclined flow field. Eq. 
(26) then is generally written as 
der + (7/17) + (Hew/4?) — (M? — 1)brz = 0 + 069?) (27) 
where for supersonic flow M? — 1 = B* > 0. 
If ¢(x, r) is a solution of Eq. (27), then Lamb* has shown that 
o(x, r, w) = cos w[d¢(x, r)/Or] (28) 
is also a solution of Eq. (27). Now Eq. (6) is a solution of Eqs. (4) or (27), and by comparison with Eqs. (8) and 
(28) a solution of Eq. (27) may be written 


cos he Br [fo(t) — fe(x — Br) |d =o (x/r)fo(x — Br) 
4dr A [(x — 6)? — Bet}3/2  /y2 — By \ 





o(x, 7,0) = — (29) 


This is the desired solution corresponding to a supersonic “‘doublet’’ distribution and may be used to represent the 
effect of the transverse flow component of the inclined uniform flow. The boundary condition to be satisfied by 


the transverse flow alone is 


(x, r, w) + v.. cosw + r’¢,(x, 7, w) = 0 (30) 


The effect of the transverse flow (v. = au.) and “doublet” distribution alone may be found and then super- 
imposed on the previously found axial flow solution, since Eqs. (4) and (27) are linear. However, without addi- 
tional justification it does not seem logical to consider any Mach Number effect due to the combined flow, since 
the Mach Number for the axial flow is u.. /a while that for the transverse flow is essentially au,./a in this type of a 


superposition. 
J 


Eq. (12) shows that Eq. (29) may be written as ' 
’ 


: B? . 
(x, 7,0) = oe io + 5 ( 108") fel! (x) + 0(Brf) | (31) 
T r é - 
and therefore 
, B?2 
Ay ete his bara M 5 (10 og r) i") + oBrf) | (32) 
ole, ra) = EY — BO) 5 A liog®) fia) + 018% | (33) 


Substituting Eqs. (32) and (33) into Eq. (30) shows that 
fr(x) = trattar*| 1 + B* (1og” (rr’’ + r’*) + O(r?) | = 4au, st + B? (ios ) |+ O(r*) (34) 


Then Eq. (32) becomes 
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Aus - 2(X, 7, w) ame 2ar'| 1 + = (loe’ (=F = Ss” )| a 
Ue Ug COS w& T x Ss 2 


2,00? 
2ar'| 1 + B? (10. )(r 4+ er’? 4 7 > )| + O(r*) (35) 





and Eq. (33) becomes 


Ave/us = (x, 7, w)/U. COSw = —a + O(r? log r) (36) 

The velocity components due to the axial flow alone are given by Egs. (18) and (19) as 
Au/u,. = (r’? + rr’’) log (r/x) + O(r?) (37) 
Av/u.. = r’ + O(r? log r) (38) 


The static pressure is given by Eq. (26) as 


a2 7=2 ae 1 my Pope 
bo (Sy =f al etrearner- tal 
Po Le 2a.” -- 


as A A . Av Ar F 
{i+ Starlite (1 + 4 cosa) — (acoso + 4 cos a) — 
2 u 


dbu\? y/(y¥-1) 
(- a sin w +%) } (39) 
r 
Then 


2 (#1) = —{fo(#)+(#) +a — a(t) “ J( =e 
CG = (2 1) 2 Uo " u sl Cl ai Ht 2 - . 7 ” Un ie 
(«+ eye fos w + [ca — a( TY + (JE) + 203 fost — 2ein of + 
Us /Ua Ue “ 


O(a?, Au*®, Avt) = 4ar’ cos w — 2(Au/u,)—(r’)? + O(rs,a7) (40) 








The resultant transverse force per unit axial length is given by 


. Av2\ A 
F, = —- i gC,r cos w dw = 2xrq (2 + = rae (« + an) a (41) 
0 Uo Uo Uo 


Then, substituting Eqs. (35) to (38) into Eq. (41) yields 
F, A x B? o,f 
am oer (i+< + =X ~) + O(r*) = 4rarr’ {! + log (: iG 4 1)r’? + (4B? + Aer! + alt (42) 
Un u 


The lift and moment can then be computed from 


C, = (1/gA) Sc” F, dx t (43) 
= (1/gAL) fo” F,x dx| 


Eq. (42) shows that if the cross product term Au X Au: is neglected then the lift and moment must be obtained 
only by 


F,/q = 4rarr’ + O(r‘ log r) = 2aS’ + O(r! log r) (44) 


which is the expression originally obtained by Munk’® for the incompressible potential flow about a very slender 
body of revolution. This also proves that the expressions given by Tsien? are incorrect, since they entirely neglect 
Au X Au, while extending the order of approximation of Au alone. 

Lighthill* has proposed that only Eq. (44) is permissible in these linearized supersonic flow solutions. It will 
now be shown that this is true, at least for the flow about an infinite cone. 

Ferrari’® has shown that, for an infinite cone of semivertex angle 0, the solution of Eqs. (11) and (30) is 


Bau. cos w [xv (x? / Br?) - =-1 — Br cosh~ 1 (x/ Br) 


(x, ’, w) = — jap 


= : (45) 
(B/tan eV (1) 1, B? tan? 8) — 1+ B?* cosh (1 /B tan 6) +2V1 — B* tan? 0 





where r = 6x. If, instead of Eq. (30), the first order boundary condition 


o,(x, 7, #) = —au. COS w (46) 











638 JOURNAL OF THE AERONAUTICAL SCIENCES—NOVEMBER, 1947 








is used, the resulting expression for the potential is in agreement with Eq. (45) only if the order of the solution is 
restricted to 6 = 0. To this order Eq. (45) becomes 


(x, 7, w) = (au. cos w)x6(1 + 26°B* log 6) + mf 


¢:(x, 7, w) = (2au.. cos w)6(1 + 6?B? log 6) + 0(6*) (45a) 


—¢,(x, r,w) = au. cos w + 0(6?) 


However, these expressions are the same as given by Eqs. (31) to (36), indicating that, in general, if the first order 
boundary condition Eq. (46) has been used then only Eqs. (35), (36), and (42) should be used. If Eq. (45) is sub- 
stituted into Eq. (26), however, it is satisfied only to the order 0? = 0, showing that, at least for the infinite cone, 
only Eq. (44) is valid for the transverse force distribution. Asa matter of fact, Eq. (44) should be the only rigorous 
linearized solution of this type for any body, since it is not logical to superimpose an axial and a transverse flow 
unless the Mach Number effect is analogous in both cases. 

If the reference area A is the cross-sectional area of the base of the sharp nose body of revolution, then Eqs. 
(43) and (44) yield the following simple relations: 


L 
dC, —. dn _ 2(1 ¥ e =“) ‘ 2(1 ” a) (47) 





° eae * ae AL AL 
Fig. 4 compares Eq. (47) with the experimental data from Géttingen'! for two bodies of revolution at an angle of 


attack of 3°. 
The pitching moment and resultant lateral force due to a pitching velocity may also be obtained in a similar 
manner. Ifthe pitching velocity is °, then Eqs. (32) and (33) may again be used with the new boundary condition 


(x, 7, w) + & (x — Xq-) cosw + 1’g,(x, 7, w) = 0 (48) 
This results in a transverse force distribution of 


4 = — [4arr'(xae — xX) — Qari] + O(r4 log r) = —2[(xee — x)S’ — S] oo + O0(r‘ log r) (49) 


° 


@ 


which agrees with that given by Munk® for incompressible flow. 

In a circular flight path or with an instantaneous pitching velocity (a), the body of revolution with a reference 
base area of A will experience, in addition to the lift and moment given by Eq. (47), the additional lift and moment 
given by Eqs. (43) and (49) as 
aC, 2L 
= — [1 yee (xac/L)] 


2 f° (MMOH OOD: " 


2L ee moment of voleme tad vue | 
Use AL? - AL } 


Xee/L = ¥Cy/> Cz. 


The boundary condition, Eq. (48), is dependent on the location of the aerodynamic center, Xa, from Eq. (50), since 
the resultant transverse force vector must pass through the instantaneous center about which the body is rotating 
and the moments acting on the body must be in equilibrium about the point in the body where this resultant force 
vector intersects the body axis. 

Eqs. (49) and (50) are independent of the Mach Number. The Mach Number would not enter these equations 
unless terms of the order r‘ log r are considered in combining Eqs. (32), (33), and (48). However, by analogy with 
the development of Eq. (44), it does not seem possible to obtain any Mach Number effect on the transverse force 
from the linearized theory. The higher order additional terms, which are implicitly contained in Tsien’s solution, 
contain B* and therefore excessively decrease the lift and pitching moment coefficients with an increase in Mach 
Number beyond +/2. However, by considering the development of Eq. (20a), it is seem that Eqs. (44), (47), 
(49), and (50) should only be applied for M = +/2. 








where 


THE SUBSONIC AXIAL FLow ABOUT A Bopy OF REVOLUTION 


Proceeding in an analogous manner, the author has shown in reference 12 that the axial subsonic flow about a 
sharp-nose slender body of revolution is represented, for the same body shape in either subsonic or incompressible 
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flow, by 


ue  4rJo [(x — 8? + BY)?” 
Av B% ft S'(#)dé 


uw 4rJo [(x — 8? + 6%] 





+ O(r* log r) 


Au ef" S'(E)(E — x)dé 


(51) 








a. O(r* log r) 


where 8 = 1/1 — M? S11. 

Eq. (51) may be used as a first order approximation for the subsonic flow on or about a slender sharp nose body 
of revolution or as a solution, exact to the order r‘ log r = 0, for the incompressible flow (6 = 1) about any body of 
revolution. For example, Eq. (51) shows that the wind-tunnel wall interference correction would be given through 
the use of 


Au/u., = (1/446*r*) fo S"(t)(x — é)dé 


so that the wind-tunnel wall compressibility correction would vary with 1/(1 — M*)*’*. 
In reference 12 it is shown that Eq. (51) may also be written on the surface of the body of revolution as 


du _ _ iffi- Ce ES O's V5( *)) ay =) Mn 

a, 34 ee 2x (1 + tog F ee aN a +3(1 ep + | 
2 9 Dy? " 

: (1 _- = + 2) so > abt 0(r3) (52) 


24 | oe * 


/ 


Av/u.o = r’ + O(r? log r) 


where the surface pressure coefficient is again given as in Eq. (20) by 





7 ) 
1 ae Br L 2x 
=({- pS —- ” =f ij vy? 
C, dl : x Ox 1 + log am 42 (1-2) 4 SS” + <( L )s + | 
ae L } 4 ” 
+ (53) 

L? Qe 2x? 
sil! ms r+ )s +--+. [= (1)? ) + O(r*) | 
“ - 


It is obvious that for the case of subsonic flow all of the terms in Eq. (53) may be justified as an approximation for 
the surface pressures on a body of revolution in compressible flow. This is true because, when 6 = 1, Eq. (53) 
is a solution of the exact potential equation for incompressible flow to the order r* = 0. For 8 < 1, the Mach 
Number effect is included as log 8 in the terms of order r* and therefore may be retained for 8 > 0 to give a valid 
approximation for the effect of Mach Number. In the supersonic case the terms of order r? contain log B, which 
can only be retained for B values not much greater than unity, otherwise these terms will become numerically 
greater than the lower order terms as the Mach Number increases. 

Eq. (53) provides the surface pressure distribution on a slender symmetrical body of revolution in either in- 
compressible or subsonic potential flow. The first order compressibility effect is only valid for a slender sharp- 
nose body of revolution, but for incompressible potential flow Eq. (53) is a valid approximation even for a rounded 
nose symmetrical body of revolution, as shown in Fig. 5 where Eq. (53) is compared with the exact solution for a 
prolate spheroid of 10 per cent thickness ratio. 

Eq. (53) shows that, unlike the simple first order compressibility correction of 1/8 for two-dimensional subsonic 
flow, the first order Mach Number correction for a slender body of revolution varies with log 8, and the magnitude 
of the correction is affected by the body shape, as well as the actual location on a given body. Fig. 5 shows that 
for a prolate spheroid the correction varies only slightly with the location on the body until the stagnation points 
are approached. Of course, the compressibility correction given in Eq. (53) is not applicable near the stagnation 
point or for large values of r’; however, it may be applied as a general first order ratio correction. For example, 
the compressibility correction for the maximum negative surface pressure coefficient for any symmetrical body of 
revolution would be obtained from Eq. (53) by taking x = L/2 so that 


CG, ) # /| ae (es )] 
(Gr), = 1+ Vow e/|1 +1985 - (Gssr +--+) If (54) 


Fig. 6 compares Eq. (53) with the incompressible potential flow solution of a 10 per cent thickness ratio body of 
revolution generated by a parabolic arc. Also shown in Fig. 6 are the supersonic surface pressure coefficients ob- 
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tained from Eqs. (20) and (20a). 
distribution. 
aft portion than would Eq. (20). 


1947 


The effect of the wave drag is seen as the unsymmetrical shift in the pressure 
Eq. (20a) would, of course, predict less total wave drag and predict more flow compression at the 
The values given by Eq. (20) for this type of a body are in close agreement with 


the results of von Karman" ° and Lighthill‘ for Mach Numbers near +/2. 


THE SUBSONIC FLOw ABOUT AN INCLINED Bopy OF REVOLUTION 


For the inclined subsonic flow the doublet distribution analogous to Eq. (29) is easily found from Eq. (28) and 


reference 12 as 


(x, r,#) = —(Br cos w/4m) So” falt)dé/[(x — £)? + B*r?]}*”? 


Following reference 12, this may be written as 


Br cos : 
Ar Br 





o(x, 7, w) = 


Br log (2) | + 0(Brf) (55) 


The boundary conditions are again given by Eq. (30), and, combining this with Eq. (55), 


f2(x) 


Au 


Uo 
Ave/ tn 


—a + 0(r? log r) 


2au., BS[1 — B? log (r/L)(S’’/2x)] + O(r*) 
Qtr 
2ar’ E — B* log (F)(r" + 49? + * - )| + 0(r*) (56) 





From Eq. (52) the first order terms for the axial subsonic flow components are given by 


Au/ue = (S'’/2r) log (r/L) + Ts (57) 
Av/u. = r’ + O(r* log r) ‘ 


Substituting Eqs. (56) and (57) into Eq. (41) yields a resultant transverse force distribution of 


q 


The first order approximation is again given by the 
simple Eq. (44), showing that Munk’s® expression for 
the transverse force distribution is a valid first order 
approximation for supersonic, subsonic, or incompres- 
sible flow. Eq. (58) is valid only for sharp nose bodies 
of revolution of such small fineness ratio that Eq. (57) 
gives a suitable approximation to the axial flow surface 
pressures. Eq. (52) is valid even for a rounded nose 
body of revolution, but the same cannot be said of 
Eqs. (57) and (58), which give erroneous results for 
these cases. However, it may easily be shown that the 
exact expression for the incompressible potential flow 
about a prolate spheroid is given by 


F,/q = 4mra rr'/(1 + r’?) (59) 


In reference 13 it is shown that Eq. (59) is in excellent 
agreement with the transverse force distribution on 
arbitrary rounded nose bodies of revolution even at 
rather large angles of attack. 

Fig. 7 compares the transverse force distributions 
given by Eqs. (44), (58), and (59) with the experimental 
values, also used in reference 13, for a typical dirigible 
shape at an angle of attack of 10° and a @ value of ap- 
proximately unity. It is seen that Eq. (59) gives a 
good approximation for the transverse force distribu- 
tion as long as there is potential flow with no separa- 
tion. Eq. (59) should be used for all subsonic flow 


2y2yttt 
Os 4rarr’ {1 — log (io — 1)r’? + (48? — 1)rr”’ + Pa -|\ (58) 


r 





computations, and perhaps the effects of compressi- 
bility could be estimated by using Eq. (58) as a ratio 
correction, although Eq. (58) is rigorously valid only 
for incompressible flow (8 = 1) since the axial and 
transverse flows have different Mach Numbers for 
8 <1. Ima real viscous fluid the flow will separate 
and not follow a closed trailing-edge contour; there- 
fore, if the wake shape is known, a much better approxi- 
mation may be obtained by using this wake shape as 
the aft portion of the body contour. If Eq. (59) is 
applied to the wake shape shown in Fig. 7, it will give 
values that are in close agreement with the experimen- 
tal data. 

For supersonic flow, however, as previously men- 
tioned, only Eq. (44) should be used as a rigorously 
valid approximation. Eq. (44) shows that a body of 
revolution with a closed trailing edge (A = 0) will have 
no lift but still have an unstable pitching moment. Of 
course, the actual wake existing in a real viscous fluid 
will tend to produce a lift increment. If the cross- 
sectional area of the wake is taken as A, Eq. (47) should 
give results agreeing with the experimental data. 

It is easily shown that Eqs. (49) and (50) also apply 
for a slender sharp nose body of revolution having a 
pitching velocity in subsonic flow, at least for 8 values 
considerably greater than zero. Munk’ has shown 
how Eq. (49) may be modified by the apparent mass 
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moment of inertia in order to approximate the trans- 
verse force distribution for a rounded nose body of 
revolution. However, for a body of revolution that is 
closed—that is, rounded or pointed at both ends—the 
resulting lift and moment of Eq. (49) will always vanish 
so that all the terms of Eq. (50) are zero for this case. 
However, the actual wake existing in a real viscous fluid 
will again modify these results, and an agreement with 
the experimental data should be made possible by tak- 
ing A as the cross-sectional area of the wake at the 
trailing edge. 


DISCUSSION 


It is interesting to note that only in the two-dimen- 
sional supersonic flow solution is the wake behind the 
closed body eliminated; that is, in a theoretical non- 
viscous, compressible fluid the flow behind the trailing 
edge of a two-dimensional airfoil immediately resumes 
the tree stream direction and static pressure (if it is 
assumed that the shock waves are so weak that the 
change in entropy may be neglected, as would be true 
for a slender sharp nose airfoil at small angles of at- 
tack’). Jones!‘ has shown that a body of revolution in 


supersonic flow has a wake similar to that which it has 
in subsonic flow. This may be shown by writing Eqs. 
(16) and (51) forr = Oas 


se) ate! Sapte vith ie (for M>1, x>L) 
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This shows that the supersonic wake is identical to the 
subsonic wake and has exactly twice the pressure dif- 
ference, and both the subsonic and supersonic wake 
pressure coefficients are independent of the Mach 
Number to the first order. The integration of Eq. 
(60) along the axis either ahead of, or behind, the body 
is, of course, rigorously valid for the case of incom- 
pressible flow (8 = 1). By analytical continuity the 
subsonic case should also be valid, at least for M << 1. 
Although there may be some doubt as to the validity 
of applying this procedure to supersonic flow, still 
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from physical considerations the integral of Eq. (60) 
seems to be valid. The two-dimensional supersonic 
flow, of course, has no wake if the flow is isentropic 
throughout, because, as shown in reference 7, the pres- 
sure depends only on the angle between the slope of the 
local velocity vector relative to the free stream direc- 
tion. 

It should also be noted that theoretically an ordinary 
stagnation point occurs at a sharp trailing edge in each 
case except that of the two-dimensional supersonic 
flow. Of course, the equations presented give an in- 
finite negative velocity increment at the trailing edge, 
but this only occurs because the points (0, 0) and (Z, 
0) are singular points that should only be approached, 
never actually reached. For example, if a symmetrical 
parabolic arc is taken as the body profile, it will be 
found that, in the two-dimensional subsonic flow with 
8 = 1, the equation in reference 3 or 15 for the source 
distribution representing the parabolic are gives a 
stagnation point at x/c ~ 1 — 10~‘; in the axially 
symmetric subsonic flow with 8 ~ 1, Eqs. (51) or (53) 
show that a stagnation point occurs at x/Z ~ 1 — 10—"}; 
and in the axially symmetric supersonic flow with B = 
1, Eqs. (16) or (20) show that a stagnation point occurs 
at x/L ~ 1 — 10~. 


CONCLUSIONS 


The surface pressures on a symmetrical body of revo- 
lution in axial subsonic flow are satisfactorily given by 
Eq. (53). For incompressible flow, a higher order ap- 
proximation is given by Eq. (51) with 6 = 1. For 
supersonic flow at Mach Numbers near 1/2, Eq. (20) 
provides a satisfactory solution for uniformly continu- 
ous bodies of revolution other than those having a coni- 
cal shape, for which Eq. (20a) should be used. The 
Karman-Moore! solution may be used for bodies of 
revolution that have discontinuities in the derivatives 
of the cross-sectional area, but the computations should 
be made only for M ~ +/2. None of the linearized 
supersonic solutions are valid, however, for Mach Num- 
bers less than the minimum Mach Number required 
for the existence of an attached conical shock wave or 
for predicting any variation with Mach Number for 
M >> V2. 

The lift and moment acting on an inclined body of 
revolution in supersonic flow are given directly by Eq. 
(47), while for subsonic flow Eq. (59) should be inte- 
grated to determine the lift and moment. The damping 


moment due to the pitching velocity of a body of revolu- 
tion in either subsonic or supersonic flow is given by 
Eq. (50). If the body of revolution is closed at the 
tail end (A = 0), it would be best to assume a probable 
shape for the wake and separated flow region and then 
use the outer streamlines of this region (A > 0) in ap- 
plying Eqs. (47), (50), or (59). The lift and moment 
results are independent of the Mach Number, but in 
the supersonic case these results do not apply for Mach 
Numbers less than +/2, while in the subsonic case they 
do not apply as the Mach Number approaches unity. 
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The Dynamics of a Swept Wing 


H. R. LAWRENCE* 
Northrop Aircraft, Inc. 


SUMMARY 


It is assumed that a swept wing may be represented dynami- 
cally by a plane curved rod of variable cross section with finite 
stiffness in torsion and in bending normal toits plane. The equa- 
tions of motion for this configuration are exhibited with the re- 
striction that the curvature of the elastic axis be finite. By 
allowing the curvature to become infinite at a point, certain 
auxiliary conditions are obtained at the bend. The validity of 
this limiting process is examined from the physical point of view. 
A wing with an elastic axis consisting of a plane broken line of 
(n + 1) segments is denoted as an ‘‘n bend wing.”” A method is 
developed for calculating the natural frequencies and modes of a 
single bend wing. Formulas for the natural frequencies and 
modes of a free-free symmetric single bend uniform rod are ob- 
tained. A generalized type of orthogonality condition between 
the modes of a single bend wing is proved. Based on these re- 
sults, formulas are developed for the dynamic response of a single 
bend wing to a general system of forces. A method of generaliz- 
ing the above results for a single bend wing to the m bend wing is 
indicated. The use of an m bend wing as an approximate repre- 
sentation of a wing of continuous curvature is discussed. 


NOTATION 


= length of rod measured along elastic axis 

= bending stiffness of wing 

= torsional stiffness of wing 

= r-coordinate of section center of gravity 

external running shear load 

= external concentrated shear load applied at s = a 

= external running torque 

= external concentrated torque applied at s = a 

= mass moment of inertia per unit span about section 
center of gravity 

= same about section elastic axis 

= mass per unit span 

internal bending moment 

= external concentrated bending moment applied at 
s=a 

= see Eq. (16) 

= distance measured normal to elastic axis 

distance measured along elastic axis 

time variable 

internal shear 

internal torque 

kinetic energy 

vertical displacement at any point 

= vertical displacement at the elastic axis 

= potential energy 

x, y,2 = Cartesian coordinates 
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a = angle of sweep 
8 = angle of twist 
m = (B/C) at a bend in the elastic axis 


The subscripts m, n denote integers in the following symbols: 
Amn = Bun = see Eq. (32) 
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C, = see Eq. (29) 

H, = see Eq. (39) 

k,* = mw, */B for uniform beam 

1,2 = Jw,?/C for uniform beam 

P, = bending displacement of uniform beam in mth mode 

Qn, = angle of twist of uniform beam in nth mode 

W, = natural frequency of mth mode of uniform beam 

¢n = bending displacement of nonuniform beam in mth 

mode 

vn = angle of twist of nonuniform beam in mth mode 

An = natural frequency of nth mode of nonuniform beam 

T'mn = see Eqs. (35) and (41) 

on = see Eq. (48) 

En = see Eq. (47) 

Note: The subscripts s and ¢ are used to denote partial dif- 
ferentiation. 


INTRODUCTION 


bbw USE OF WINGS with large sweep has in the past 
been restricted to airplanes of the flying-wing type. 
With the advent of airplanes designed to fly in the 
transonic and supersonic régimes, the highly swept 
wing appears destined to become of increasing impor- 
tance. 

The problem of calculating the dynamic response of 
an airplane with a swept wing may be considered in 
three parts. First, a simplified model of the structure 
is chosen. Second, a method of calculating the dy- 
namic response of the simplified model is developed. 
Finally, evidence is assembled to show that the sim- 
plified model represents the airplane dynamically. 

In the case of an unswept wing, the two most fre- 
quently used models are: (1) a set of spring-connected 
masses and (2) a rod of variable cross section. 

A generalization of the second dynamic model noted 
above to a swept wing is the subject of this paper. 

It is assumed that a swept wing may be represented 
dynamically by a rod of variable density and cross sec- 
tion possessing an elastic axis in the form of a plane 
curve. Only loads and moments normal to the plane 
containing the elastic axis are considered. The fuse- 
lage and nacelles are considered as affecting the mass 
and stiffness characteristics of the wing without intro- 
ducing additional degrees of freedom. 


I. GENERAL EQUATIONS OF MOTION 


It is convenient to consider first a rod whose greatest 
diameter of cross section is small compared with the 
radius of curvature of its elastic axis. It is assumed 
that the dimensions of the wing normal to a plane con- 
taining its elastic axis are small so that the wing may be 
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Fic. 1b (bottom). 


Fic. la (top). 


treated geometrically as a plane figure containing an 
elastic axis. The elastic axis BAPB’ is shown in Fig. 
la lying in the plane OXZ. 

It is evident from the figure that any point on the 
undistorted wing may be specified uniquely in terms of 
the coordinates s and r. Let u represent the displace- 

— 


If v represents 
ed 


the displacement of the point P in the direction OY 
and £ represents the rotation of the line PQ about the 
elastic axis at the point P, then u = v + rf. 

Consider a typical slice of the wing taken normal to 
the elastic axis as illustrated in Fig. lb. The line PP’ 
is the elastic axis where BP < BP’. The external 
loads are represented by a load per unit span f and a 
torque per unit span g. The internal elastic forces are 
represented by a shear S, torque 7, and bending mo- 
ment M. 

The equations expressing the internal elastic moments 
in terms of the deformations of the structure may be 
written! as follows: 


ment of a point Q in the OY direction. 


M = Biv, — B/p), 
T = C(B, + 2,/p) (1) 
> = —M, 


where B and C represent the bending and torsional 
stiffness, p represents the radius of curvature of the 
elastic axis in its plane, and subscripts are employed 
to denote partial differentiation. 

The equations of motion will be deduced from Hamil- 
ton’s Variational Principle. By using this procedure 
the proper end conditions are obtained in a purely 
mechanical fashion. The kinetic energy 7 and the 
potential energy V must be expressed in terms of the 
mass, stiffness, and geometric parameters. 

It is evident that the kinetic energy of a slice of the 
beam, made by planes normal to the elastic axis and a 
distance ds apart measured along the elastic axis, is 
given by the expression 


dT = '/2[m(v, + eB,)? + JoB7]ds 
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where m is the mass per unit span, e is the r-coordinate 
of the center of gravity of the slice, and Jo represents 
the mass moment of inertia per unit span of the slice 
about its center of gravity. Upon integrating over the 
span and defining J as the mass moment of inertia per 
unit span about the elastic axis, the kinetic energy may 
be written as 


T ='/2f," [mv? + 2mev, 8B, + JBZ\ds (2) 

The potential energy V may be written as the sum of 
two parts. Let Vi represent the elastic energy and JV, 
represent the potential energy of the external forces, 


The internal potential energy may be written in the 
familiar form of 


Vi = "/2 ft" [((M?/B) + (T?/C) ds (3) 


where M and T are expressed in terms of displacements 
by means of Eq. (1). 
The potential energy of the external forces may be 
expressed in the form 
V2 = —Si' (of + Bg)ds (4) 
Hamilton’s principle states that the integral A be- 
tween arbitrary fixed limits /, 4, where 
A= Jf," (T — V)dt 


is an extremum for the motion or 6A = 0. 
Upon performing the variational process and collect- 
ing terms, the following equations are obtained: 


MV i, + meB 1 = —M,, + (T/p)s +f, l ae A (5) 
mev;, + JBu = T; + (M/p) + g, ftzeren 


Miv, = 0 
[S + (T/p)]év = 07 at each end of (6) 
[T — (M/p)iB = 0 the beam 


II. THE EFFECTS OF A BEND 


The equations of motion obtained in Part I are valid 
only when the radius of curvature of the elastic axis is 
large compared with the largest diameter of the beam 
cross section. The effects of a vanishing radius of 
curvature are investigated first from the purely formal 
point of view. The general validity of the results so 
obtained will be considered later. A bend in the rod 
axis is considered as the limiting case of a circular arc. 
Fig. 2 represents the geometry of the elastic axis in the 
neighborhood of a point of large curvature. 

The segments B’A’ and AB are straight lines whereas 
the segment A’OA’ represents a circular are of radius 
a. The equations of motion valid over B’OB may be 
written from Eq. (5) as follows: 


mv, + meB,, = S, + (T/p)s + f\ (7) 
mev;, + JByn = T,; + (M/p) + gf 


The first equation is multiplied by ds and integrated 
from B’ to B;_ thus, 


[S + (T/e)] |B = Sa? (mou + meBu — fds 
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or noting that p = ~ at Band B’ and defining S, and Sp 
as the shear at B and B’, respectively, 


Si — So => Sz (mv os meB i — fds 


If the points B — A and B’ — A’ while a — 0, the length 
of the path of integration approaches zero. Since the 
integrand is assumed to be finite 


Si = So (8) 


Consider the equations of motion on the arc A’OA. 
The integration of the first of Eqs. (7) may be written 


as 
-M, + (T/p) = So® (moun + meBu — f)ds + const. 


or in polar coordinates 


-Mo + T = a2 fo'(moy + meBu — f)dd + 
(a) X (const.) 


Since both the constant and the integrand are assumed 
to be finite, the limit of the above equation as a — 0 
is 


In the same way the second equation may be written 
as a — OQ in the form: 


™%+M=0 


if (mev + JB) is finite. 

Let Mo, To, and M,, T; represent, respectively, the 
values of M and T at@ = —aand@= a. Eggs. (9) 
and (10) may be integrated to obtain M and T in terms 
of two arbitrary constants; thus, 


= —A,sin 6 + A2 cos 0 
Mz= A;cos@é+ Aegsin 0 


Upon substituting the values of T and M at@ = —a 
and @ = a in Eqs. (11), four equations are obtained. 
Upon eliminating the two arbitrary constants A; and 
A», the following boundary conditions at a bend are 
found: 


(Ti + To) sin a — (M, c= Mo) cosa = of (12) 
(T, = To) cos a ote (M, oe Mp) sin a = 0 


Eqs. (1) may be used to obtain two additional bound- 
ary conditions at a bend as follows. On the arc A’OA, 
Eqs. (1) may be written in polar coordinates as follows: 


—(dv,/d0) + B = — 
(dB/d0) + v, = a(T/C) 


Asa — 0, the right-hand members of the above equa- 
tions approach zero, since M/B and T/C are assumed to 
be finite. Let v, and 6 be, respectively, v;,, Bo and 2,,, 
fiat @ = —aand@ = +a. 

Then Eqs. (13), with the right-hand sides equal to 
zero and the special values at 9 = =a, are identical 
with Eqs. (9) and (10) under the substitutions 1, — M 
and 8— 7. Hence, the following boundary conditions 


are obtained: 


(10) 


(11) 


(13) 














Fic. 2 


(8: + Bo) sin a — (¥, — U4) cos a = it (14) 
(8:1 — Bo) cos a + (Um, + %.) sina = 0 


The second Eqs. (13) may be transformed into the 
form: 
dv/d@ = a?(T/c) — a(dB/dé) 


Upon integrating with respect to @, allowing a to ap- 
proach zero, and defining v; and as the values of v at 
= —6 and 8, respectively, there results 


v1 = UV (15) 
providing T/C and 6, — ( are not infinite. 

Thus, the Eqs. (8), (12), (14), and (15) give six condi- 
tions connecting the loads and deformations on the 
opposite sides of a bend in the elastic axis. The equa- 
tions fall into two sets of three. One set is independent 
of the elastic properties of the system. This set, Eqs. 
(14) and (15), may be interpreted simply as conditions 
of continuity of the beam elastic axis. The second set, 
Eqs. (8) and (12), represent the equilibrium between 
the shears and moments acting on the sides of the bend. 
In order to evaluate Eq. (12) in terms of the deforma- 
tions of the beam, it is necessary to know the ratio 
B/C. = watthe bend. The theory of elasticity as em- 
ployed herein is not adequate to produce this ratio, 
since its assumptions have been violated. On the other 
hand, it is evident physically that this ratio exists, 
since one can imagine an experiment that would yield 
it. Although the author is unable to prove that the 
only effect of allowing an infinite curvature of the 
elastic axis is to produce an undetermined elastic 
constant, results of calculations seem to support such a 
view. 

It is evident that any wing elastic axis may be ap- 
proximated by a set of straight-line segments meeting 
at bends. The simplification results from the fact 
that a differential equation with variable coefficients 
such as Egs. (5) is replaced by an equation with con- 
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The complicated coupling between bending 


and torsion exhibited in the boundary conditions, Eqs. (12) and (14), indicates that the use of uncoupled 


modes for a highly swept wing may not prove satisfactory. 
It is convenient to have a compact notation for the boundary conditions at a bend. 


fined in the following manner: 


O (x, B, s, h) 28 [ o(s + 0) 7 v(s 9) 


Uess(S + 0) aad Vess(S “si 0) 





The operator O is de- 


[B(s + 0) + B(s — o)] sin a — [v,(s + 0) — v,(s — 0)] cosa (16 
[B(s + 0) — B(s — 0)] cos a + [v,(s + 0) + 2,(s — 0)] sina ) 
[8.(s + 0) + B,(s — 0)] sin a — [u(s)][0ss(s + 0) — v45(s — 0)] cos a 
| [8,(s + 0) — B,(s — 0)] cos & + [u(s)]fou(s + 0) + va(s — 0)] sin a] 





where the expression on the right is a matrix of six rows and one column. 


III. THE NATURAL FREQUENCIES AND MODES OF A 
UNIFORM SWEPT RopD 


An analytic solution is obtained for the natural fre- 
quencies and modes of a uniform free-free rod with a 
bend at the centerline. This solution exhibits the form 
of the results obtained by use of the bend boundary 
conditions. It is also useful in the calculation of the 
modes and natural frequencies of a free-free swept 
wing. 

The geometry of the configuration under considera- 
tion is illustrated in Fig. 3. The differential equations 
and boundary conditions may be written 


mv, + Brssss = 0, (—b i b) 
JBiu — CBs, = 0, (—b<s <b) 
Uss = VUsss = Bs = Qats = +) 


O(v, 8, 0, uw) = 0 
Upon performing the substitutions for v and 8, 


vy = Pett) 
B = Qeiwt f 


in Eqs. (17), there results 


(17) 


(18) 


For certain discrete values of w = w,, the above equations have solutions P,, Q,. 
coupled natural frequencies of the system, whereas the functions P,,, Q, together form a coupled mode of vibration. 
The displacement u of any point on the wing in the mth mode is given in 7,s-coordinates by the expression u = 


P,(s) + rQn(s). 





















































Modes of a Uniform Beam — Genii 
with 23° Sweepback. ——-Twisting 
14 Symmetric 2d Symmetric 3 Symmetric. 
€,- 2.80 E,= 4544 E 5-5 4E 

2.36) (526) = el 5.20) 
‘ | F te - ‘ 
~ yi0e) (5.50) \ See 
[% Anti-Symmetric 2% AntiSymmetric | 3° Anti-Symmetric. 
E,=3./0 E,7394 E5* 5.40 
| dele | ah ony ie 
| (499) | aay 7 || ~~-“~{s-41) 
Fic. 3. 
Pos —kRP= 0, (—bd <s< b) 
0., + 2QO = 0, (—b<s < db) 
Py, = Py = Q, = Oats = +b (19) 


O(P, Q, 0, uw) = 0 
where kt = mw?/B and /? = Jw?/C. 


The values of w, are the 


The solution of Eqs. (19) shows that the symmetric and antisymmetric modes may be considered independently. 
With k, and /, corresponding to the frequency w,, the solutions may be written as follows: 


Symmetric Modes 
Case I.—w-, = Ww = 0. 
f agk = i Q-1 = 0 


Po =s, Qo = cot a, s>0 
P,’ = —s, Qo’ = cot a, ¢.< 0 


(20) 


(21) 


Case II.—w, # 0, then w, may be evaluated from the frequency equation: 


k, (1 — cos k,b cosh k,b) 


+ yu cot? a X I, tan],b = 0 (22) 





sin k,b cosh k,b + cos k,b sinh k,,b 
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and the mode P,, Q, is given by the formulas: 
2 (1 — cos k,b cosh k,,b)(cosh k,s + cos k,x) — sin k,b sinh k,,b(cosh RnX — cos k,x) 4 

















P, sin k,b cosh k,b + cos k,b sinh bb 
sin k,x + sinh k,* (23) 
Q, = 2k,(cot a/cos 1,b) cos 1,(b — s) 
Antisymmetric Modes 
Case I.—wy = 0. 
Po=s, Qo = —tana ee, 9 
= 5, Q = tan a s<0 (24) 
Case II—w, ~ 0. The frequency equation is: 
k,,(sin k,b cosh k,b + cos k,b sinh k,b) , 
2 ] = 95 
1 + cos k,b cosh k,b + ston? a Xf tenth = 0 (25) 
and the modes are: 
sin k,b cosh k,,b (cosh k,,s — cos k,s) — cos k,b sinh k,b (cosh k,s + cos k,s) 
P, = + 
1 + cos k,b cosh k,,d 
(1 + cos Rnb cosh k,,b)(sinh k,s + sin k,s) — sin k,b sinh k,b (sinh k,s — sin k,s) (26) 
1 + cos k,b cosh k,b 
Q, = —2k,(tan a/cos 1,n) cos 1,(b — s) 


A plot of the first three modes for particular values of the stiffness parameters is given in Fig. 4. 
It is evident that an analytic expression for the natural frequencies and modes of a uniform swept rod with any 
set of homogeneous end conditions may be obtained. 





IV. THE NATURAL FREQUENCIES AND MODES OF A (1) There exists a denumerably infinite set of eigen- 
SWEPT WING values or natural frequencies \, such that A, < A, 41 
for all n. 


A method is developed for calculating the natural 
frequencies and modes of a free-free swept wing. It 
will be noted that the same procedure is applicable to 














other end conditions with only minor modifications. FIGURE 4~COUPLED MODES OF A UNIFORM BEAM 
The configuration geometry is again illustrated by _ 
Fig. 3. The differential equation and boundary con- a GB + .228 d+ 235° 
ditions may be written as follows: _. Slee Wchanicraas in/eS b= 161 es 0 
MU, + meBy, = —(Boss)ss, (-8 <3 < ” (27) SYMMETRIC ANTI- SYMMETRIC 
mei, + JBu = (CB;)s (—b<s <b) , - 
LA nel 
(Buss)s _ (Buss) — CB, = Qats= +b ° > eee anges ae 
O(, B, 0, uw) = 0 er : " 
or with 
v = e* o(s) m2 ns2 
— 
B = e# y(s) ed a 
1 = 








(Béss)s = (mp + me), (—b<s <b) 
(Chs)s = —M(mep + Jp), (—b<s <5) 


9 
(Bgss)s = Boss = CB; =@Qats = +b (28) 





o*3 
0 (¢, ¥, 0, u) = 0 > oh 
Eqs. (28) are of the generalized Sturm-Liouville : ve s —<$—— os 
type. The following results concerning such systems \ is VA rect 

may be borrowed from the mathematical litera- 
ture :? 
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(2) Corresponding to each value X, there exists an 
eigenfunction or mode (¢_, p). 

(3) The eigenfunctions or modes form a complete 
set of functions. 

Upon applying the above theorems to Eqs. (19) and 
noting that (P,, Q,) and (¢, y) satisfy the same bound- 
ary conditions, the following result may be stated. The 
function gy and y may be represented in the form: 


N 
(8) =D CaPa(s) 


N° (29) 
¥(s) = 2d CiQn(s) 
where the C, are independent of s. By choosing N suf- 


ficiently large, the above representations of ¢ and y 
may be made as accurate as required, providing all the 
modes of frequency less than Ay are used. 

Although the functions (P,, Q,) are seen to be a satis- 
factory means of representing (y, y) from the mathe- 
matical standpoint, they are by no means unique in this 
respect. An additional argument for their use lies in 
the fact that (P,, Q,) and (¢, ) satisfy twelve boundary 
conditions in common. 

The differential equations of Eqs. (28) may be shown 
to be equivalent to the problem of determining an ex- 
tremum of the expression 


L= 1/s (Boas? + Cy;”)ds = 
1/2 J" ,(me? + 2megy + Jy*) ds (30) 


Upon substituting the expressions for g and y from 
Eqs. (29) and (30), the following result is obtained: 


, ee 
L= 1/, > > (Amn say \?Bmn) CuCm (31) 


n=0 m=0 


where 
te J [>Ge) ae) * | 
CE)B)] +} 


Buy at J? [mPnPn + me(PrQn + PQm) + 
JQmQn|ds 


It will be noted that Amn = Anm and Bry = Bam. In 
order to make LZ an extremum, its derivative with re- 
spect to C, is equated to zero or 


N 
oL 
—a Amn — *Bnn)C, = 0, 
aC, > ( )G, = 0 


n= 0 
m=0,1,2...N (33) 


The above set of homogeneous linear equations will 
have solutions other than C, = 0 for all , only for 
\ = \, where 4, is a natural frequency of the wing. 
The corresponding values of C, may be substituted in 
Eq. (29) to obtain the r‘ mode (¢,, ¥;). 

If the wing is symmetrical about its centerline, the 
Eqs. (33) break up into a symmetric and an anti- 


symmetric set. In each case the range of integration 
of Eqs. (32) may be taken from zero to b. In consider- 
ing the symmetric (antisymmetric) modes of the wing, 
only the symmetric (antisymmetric) modes of the uni- 
form beam need be employed. 

In the case of the antisymmetric modes, the A,, 
and Bn, are calculated from the modes given by Eqs, 
(24) and (26). Since the terms A,,, = 0, the Eqs. (33) 
take the form: 


N 
=11( Bue + > BuCx) =0 
n==1 
N 
—*BrnoCo + >> (Amn — 2 Bmn)C, = 0 
n=1 
m=1,2,...N 


(34) 


If \ = O, there is a solution corresponding to the 
rigid body motion Cy = 1C, = 0," > 1. 

If X ¥ 0, Cy may be eliminated from Eqs. (34) and 
the equations may be written as 


N 
C= — u BonCn/Boo (35) 
N 
Ze (Amn ~~} Tn) Cy = 0; m=1,2,...N (36) 
n=1 
- ! 
1 [Boo Bno (35.1) 


Tun = Bo |Bo Bra 


Since the matrices [A,,,] and [I;,,] are nonsingular, 
Eqs. (36) may be solved by matrix iteration methods 
and Cy may be determined by Eq. (35). 

In the symmetric case, the A», and By, are calcu- 
lated from Eqs. (20), (21), and (23). Since Eqs. (20) 
and (21) represent rigid body motions with zero fre- 
quency, A-1,m = 0 and Ao, = O for all m. Hence, 
the Eqs. (33) acquire the form: 


N 
-(B-, —1C-; + B-1, 0Co + u B-1, Cr) = 0 (37) 


N : 
—1( By —1C-1 + Bo, 0Co + > Bo, sci) = 0 (38) 


—)? (Bn, -1C-1 + BnoCo) + > (Amn — *Bmag)C, = 0 
m=22,3... (39) 


When A = 0, there are two solutions corresponding to 
rigid body motions C-; = 1, G = C, = 0 for all n, 
and C-; = 0, Co = 1, C, = 0 for all n. 

When \ + 0, Eqs. (37) and (38) may be solved for 
C-, and Cy: 


5 
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Upon substituting the values of C) and C, from Eqs. minor modification for most of the common end condi- 


(40) in Eqs. (39), the following result is obtained: tions. The geometry is illustrated in Fig. 3. The dif- 
nw ferential equations and boundary conditions may be 
> (Ann — *T'nn)Ca = 0 written as follows: 
n=l 

mv, + meBy, = —(Brss)ss +f, (—b <5 < 8) 


B-1,-1 B-10 B-i. 

Bo, -1 Bo, 0 Bo, m 

Bo as | Bn, 0 Ba, m 
IB-1,-1 B-1, 
ia, \ Bes 


> (41) mevy + JBu = (CBs), + 8, (-—b <s <b) (42) 
(Buss)s == Buys = CB, =0 ats = =) 
O (v, B, 0, ) =0 











} 








It is assumed that the natural frequencies and modes 
of the dynamical system are known. The differential 
equations satisfied by the modes may be written in 
matrix form as follows: 


Since the matrices [A,,,] and [T;,,] are nonsingular, 
Eqs. (41) may be solved by matrix iteration and the 
results may be substituted in Eqs. (40) to find Cy) and 
C-1. 

(leet]-sesIfa} 
V. DyNamics oF A SWEPT WING — - 


A method of calculating the response of a swept Eq. (43) will be rewritten for the rth mode as follows: 


free-free wing to a system of time varying loads is de- 
veloped in terms of coupled natural frequencies and {Byres es | _ y.2] ™ mel |e (44) 
modes. The procedure may be employed with only {—C(Wn)s}s ae vy, 





Upon multiplying Eq. (43) on the left by the matrix [¢,, y,], multiplying Eq. (44) on the left by the matrix 
[¢n ¥n], and subtracting and integrating from —d to b with respect to s, there results 


Jo, Atom val fey |= tow val Reg Tf ae = am" flo valle "7 | [fas - 
i [em Vn] ee me ll fs (45) 


By integrating by parts, the left member of the above equation may be shown to vanish if the following condi- 
tions hold: 


Gr X {B(gn)ss}s = 0 ats = +b 
B(¢r)s X (¢n)ss = O ats +b >for all m and r where n # r 


Ch-(Wn)s 0 ats = 5) 
The above conditions are satisfied for most of the common end conditions. In the case of a free-free beam 
(Gn) ss - [B(¢n)ss]s = (Wn)s = 0 


ats = +b; hence, the left-hand member of Eq. (45) vanishes; thus, 


Xe? — Az?) - m nm las ="0 46 
r ue Yrs 7] me Vn ha (46) 


If \, * X,, the integral of Eq. (43) must vanish or tHe modes (¢,, ¥,) must satisfy the orthogonality condi- 


tion. 
me || o, ds = 0, (n # 1) - 
z [vr vr] vam J Ish = by (n = 1) a7) 


From the results of the Sturm-Liouville type analysis discussed previously, the deformations v and 8 may be 
expressed in the form 


N 


(3) = D oalen(s) 
"x (48) 


a on(t)¥n(s) 


B(s) 
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Upon substituting Eqs. (48) for » and 8 in Eqs. (42) and making use of matrix notation, the result may be 


> | | ie + H (49) 


written: 


oo", "][e] - 


n=0 


1947 


Multiplying Eq. (49) on the left by the matrix [¢,, y,] and making use of the orthogonality condition [Eq, 


(46)] and Eq. (43), there is obtained: 


in [tom val], eS] = rte ftom val] eS] as + ftom val Ze 


or upon making use of the orthogonality condition [Eq. (47)], 
Gn + Mon = (1/5m) Ss (en f + vng)ds 


Let 


H,(t) = (1/¢n) Jus (gn X f + Yng)ds 


Then the above equation may be written: 


Gn + r,7o, = H,(t) (51) 


If there is a small amount of damping in the system, 
it may be taken into account approximately at this 
point by replacing Eq. (51) by 


Gn + 2yAnon + Anon = H,(t) 


where y represents the fraction of critical damping. 
A solution of Eq. (52) such that o,(0) = ¢,(0) = 0 is 
given approximately for small y by the expression 


(52) 


1 r 
o,(¢) = X f H,(r)e~ *7¢-9 sin \,(¢ — r)dr (53) 
nJo 


Upon substituting the value of o, from Eq. (53) in 
Eq. (48), the dynamic response of the wing to dis- 
tributed torque and running shear may be calculated. 
In the case of a concentrated load F(t) applied at a 
point s = a on the span, the value of H/,(¢) is given by 
a limiting process as 


H, (t) is ¢n(@) F(t)/Sn 
Similarly, for a concentrated torque G(¢t) applied at s = 
a, 

A, (t) = Ya(@)G(t)/on 
The case of a concentrated moment M(t) applied at s = 
a may be obtained by taking equal and opposite forces 
+= M(t)/2e applied at a + e€ and going to the limit. 
The value of H,,(t) becomes 


H(t) = [M(t)/Sn] [den(a)/ds] 


CONCLUSION 


A method has been developed for calculating the 
natural frequencies, coupled modes of vibration, and 


(50) 





the response of an airplane to an arbitrary time de- 
pendent set of loads. This method may be used di- 
rectly to compute the dynamic loads due to landing’ 
when the external forces are known. The flutter 
characteristics of an airplane may be investigated by 
an extension of the procedure.‘ The derivatives used 
in dynamic stability of airplane may be modified for 
wing flexibility by another extension. 

The utility of the results obtained in this paper de- 
pend on the assumption that an airplane with a swept 
wing may be represented dynamically by a rod of 
varying cross section. The statement that a rod of 
variable cross section represents an airplane with a 
swept wing may be put more precisely in the following 
manner: 

(1) The kinetic energy of the airplane may be rep- 
resented by Eq. (2) by proper choice of mass parame- 
ters m, e. J. 

(2) The -train energy of the airplane may be repre- 
sented by Eq. (3) by proper choice of an elastic axis 
and stiffness coefficients B, C. 

It is the opinion of the author that both of the above 
statements are true for many airplanes in an engineer- 
ing sense. However, methods of calculating the posi- 
tion of the elastic axis and the coefficients B and C are 
not completely satisfactory. 
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Fluid Dynamic Aspects of Axial-Flow 
Compressors and Turbines 


JOHN R. WESKE* 
Rensselaer Polytechnic Institute 


SUMMARY 


Pertinent facts of the two-dimensional theory of flow through 
blade grids are used to establish relationships for the three- 
dimensional flow through axial-flow blade rows. Certain aspects 
of flow in boundary layers and wakes encountered in blade rows 
are investigated qualitatively, particularly with regard to their 
effect upon operating characteristics. A criterion of stability of 
flow is proposed and defined mathematically. 


SYMBOLS 


= diameter, ft. 

volume rate of flow, ft.* per sec. 
total head, ft. 

static head, ft. 

= rate of total energy, ft.lbs. per sec. 
absolute velocity, ft. per sec. 
lift coefficient 

chord, ft. 

pitch, ft. 

rotative speed, 1 per sec. 

static pressure, lb. per sq.ft. 
blade velocity, ft. per sec. 

= number of blades 

adiabatic exponent 

circulation, sq.ft. per sec. 
angular velocity, 1 per sec. 
mass density, lb. sec.? per ft.* 
difference between conditions in plane 2 and plane 1 
[Q/(xD?/4)u] flow coefficient 

= Ah/(u?/2g) pressure coefficient 
= h,'/h coefficient of reaction 


ll 


zx & 6 SOS eS OS RS SS A a 
ll 


1 and 2 refer to 
and T re- 


Subscripts: wu is component parallel to u; 
conditions upstream and downstream, respectively; 


fers to tip diameter. 
Superscripts ’ and ” refer to the moving and fixed rows, re- 


spectively. 


INTRODUCTION 


re OF THE CAPACITY and efficiency of tur- 
bine and axial-flow compressor stages on the basis 
of the present two-dimensional blade grid theory lead 
to extremely high values. While corresponding test 
results have been improved greatly over earlier results, 
chiefly through the application of the two-dimensional 
theory, there still exists a large discrepancy between 
experimental and theoretical values which cannot be 
explained without consideration of factors relating to 
the three-dimensional nature of the flow, inherently 
neglected in the two-dimensional theory. It is ex- 
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pected that a better understanding of these three- 
dimensional effects will improve the correlation of 
theory and test and possibly lead to an extension of the 
present limits of capacity and attainable efficiency. 
To this end an attempt is made to analyze the flow 
through axial-flow blade rows of turbines and compres- 
sors, with particular consideration of three-dimensional 
effects. 


ENERGY TRANSFER 


The problems to be considered center around the 
process of conversion of kinetic energy into shaft work, 
or vice versa. This process, although but one link in 
the chain of various states of energy, is next to the 
question of thermodynamic availability of energy the 
controlling factor influencing stage characteristics. 
It is primarily a flow process subject to the laws of 
fluid dynamics. Basic relations may be stated as 
follows: If the flow through a stage is assumed to be 
rotationally symmetrical, an element of this stage, 
bounded by two neighboring coaxial stream surfaces, 
may be considered, for which the rate of energy transfer 
is given by Euler’s equation: 


E = pQwA(rcy) (1) 


or, when introducing the circulation about the blade 
section 


T = 2rA(re,) = (T2 — T;)/z (2) 
in Eq. (1) 
E = Qp(w/2r)2T 


Eqs. (1) and (2) apply to any coaxial stream surface. 
For pure axial flow, Eq. (1) reduces to 


E = Qp(uAc,) 


The equations may be applied to an entire stage only 
if the circulation about the blade sections in all stream 
surfaces is the same. 

Eq. (2) states that the rate of energy transfer de- 
pends upon the action of a rotating system of vortices 
upon a streaming working fluid. In view of this de- 
pendence of energy transfer upon the action of vortices 
(vortex = ‘“‘turbo’’), it seems appropriate to refer to tur- 
bines and compressors collectively as ‘‘turbomachines.”’ 
The form of the expression for the rate of energy trans- 
fer, Eq. (2), is analogous to the equation defining energy 
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VELOCITY DIAGRAM AND SCHEMATIC BLADE SHAPES 
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conversion in electrical machinery. It is customary to 
express the rate of energy transfer in terms of non- 
dimensional quantities—namely, the flow coefficient 


Q 4Q 








= “= (3a 
? (eD*/4)u "Din ) 
and the pressure coefficient 
Ah Zac. "Ditte . 
y= ; ——— s——— (3b) 
u?/2g u anD 


These nondimensional quantities relate the primary 
performance quantities—namely, the volume rate of 
flow, Q, the difference of total head, Ah, and the rota- 
tive speed, n—to the dimensions of the stage, tenuously 
given by D. 
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Substitution of yg and y into Eq. (la) leads to: 
E = p(x*/8)D*n' py (4) 


More positive relations between the significant physical 
features of the stage and its operating characteristics 
can be established for a restricted family of stages 
designated ‘‘normal stages,” for which the conditions 
of flow, whether compressible or uncompressible, are 
identical at inlet and outlet as to direction and velocity 
of flow, except for a difference in head. 

For normal stages, one additional quantity is needed 
to establish the rational relations desired, since it is 
necessary to account for the percentage of variation of 
static head in the fixed and in the moving row which 
otherwise would remain arbitrary. A suitable non- 
dimensional quantity relating to this degree of freedom 
of normal stages is the degree of reaction, defined as 
the ratio of the difference of static head of the moving 
row to that of the stage, x = Ah,’/ Ah. 


RELATIONS BETWEEN NONDIMENSIONAL PARAMETERS 
AND THE VELOCITY DIAGRAM 


It is significant that the quantities y, y, and x appear 
explicitly in the velocity diagram of normal stages, 
see the lengths gu, yu/2, xu, in Fig. 1. In fact, these 
three quantities define the velocity diagram completely. 
As shown previously by Traupel,! they may be used as 
parameters to establish the direction of flow at inlet 
and outlet of both the fixed and moving rows of the 
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normal stage, hence, the ideal blade angles for any given 
operating condition of a blade element. 

A tabulation of the velocity diagrams and of the 
shape of blade of the fixed and moving rows, in sche- 
matic presentation, is given in Fig. 2 for various values 
of the pressure coefficient and the degree of reaction 
and for a constant value of the flow coefficient, since 
the latter is of secondary importance. The range of 
variation of both y and x is believed to cover the entire 
field of normal operation of ‘‘normal stages’ of both 
compressors and turbines. The blade combinations 
derived for two-dimensional flow, and shown schemati- 
cally in Fig. 2, constitute the elements for a functional 
design of blade rows. This figure will be used to illus- 
trate certain three-dimensional aspects of flow in blade 
rows. 


VARIATION OF CONDITIONS OF FLOW AND OF SHAPE OF 
BLADES ALONG THE RADIUS 


Since the blade velocity of a moving axial-flow blade 
row and also, in general, the circumferential velocity 
component of the working fluid varies with radius, 
there will be corresponding changes in the fundamental 
parameters of the stage element necessitating changes 
in the velocity diagram and changes in the shape of the 
blading along the radius. A properly designed blade 
for an axial-flow stage therefore will represent a combina- 
tion of various blade elements, each selected to meet 
precisely the functional requirements. 

Let us consider a condition of operation of a stage 
for which the axial velocity and the energy transfer are 
constant along the radius and, furthermore, for which 
the total energy upstream of the grid is constant. Then 
it follows that the total energy downstream is likewise 
constant and the flow is irrotational throughout. In 
that case the fundamental parameters vary with the 
radius as follows: 


v ¥r(rr/r)* 
Y ¢r(rr/r) 
1 — x = (1 — kp)(rr/r)? 


Fig. 3 shows the variation of the two parameters y 
and x with radius and likewise the change of shape of 
blades from tip to root in schematic presentation. Of 
interest is the pronounced change of both the pressure 
and reaction coefficient. The above applies to what is 
commonly known as a vortex type of stage. 

Let us consider further a stage of constant axial 
velocity which produces a constant increase of energy 
along the radius but for which the degree of reaction is 
arbitrarily prescribed. We choose x = '/2 constant 
along the radius. This condition is represented by Fig. 
4. It is seen that in this case the total energy / varies 
with the radius. In fact, the flow corresponding to 
conditions halfway through a row is identical with a 
solid body rotation ¢, = u/2 = const. X 1, the fluid 
rotating at half the angular velocity of the moving row. 
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It is possible to design stages for other variations of the 
coefficient of reaction along the radius which may have 
practical advantage, for instance, with regard to tip 
and root leakage of the moving and fixed blades, re- 
spectively, if the condition of radial stability is satis- 
fied, a point to which reference will be made later. 
Calculation of the variation of lift coefficient of the 
moving and fixed blade sections along the radius, 


cr! = t/ly/V o? + &) 
and sting, bastatielhe 
a” = (t/¥/Ve? + (1 — «)] 


respectively, for the conditions of Figs. 3 and 4 shows 
that, in general—e.g., for constant chord, / = const.— 
the maximum lift coefficient occurs at the root radius. 
Attempts to relate the stalling limit of the root sections 
of a blade row to a corresponding factor limiting the 
operating range of the stage have not been successful, 
and it is evident that in this respect, as in others, ac- 
count must be taken of the properties of the real fluid. 


BouNDARY-LAYER EFFECTS IN BLADE Rows 


Disregarding various other factors of importance, 
such as compressibility, attention will be given to those 
effects resulting from the viscosity of the working fluid 
which cause radial displacements of portions of the 
working fluid with respect to the main body of the fluid 
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Fic. 5. Schematic representation of flow through a curved 
channel. 


and, consequently, render invalid to some degree the 
basic assumption of the two-dimensional theory, ac- 
cording to which the flow proceeds on coaxial stream 
surfaces of revolution. 

A number of characteristics of the real fluid can be 
analyzed qualitatively, and their effect can be traced 
upon the performance of axial-flow blade rows. A 
quantitative theory for these effects, however, does not 
exist as yet. At this stage, application of the results of 
investigations of curved ducts to the flow through blade 
rows appears most profitable. It is, therefore, proposed 
to consider the flow through a blade grid as a flow 
through individual channels. 

For purposes of comparison, the incompressible flow 
through an elbow of rectangular cross section will be 
analyzed (Fig.5). It is assumed that the flow approaches 
the elbow through a straight duct of length sufficient to 
establish appreciable boundary layers. At the inlet 
of the curved channel, centrifugal acceleration produces 
a displacement of the main body of the fluid toward 
the outer wall (front of the blade) and a displacement of 
the surrounding fluid in the boundary layer, having a 
lower kinetic energy, toward the inside wall (back of the 
blade). An increase of the boundary-layer thickness 
along the back of the blade is observed. The displace- 
ments of fluid give rise to velocity components in the 
plane normal to the axis of the channel, which are di- 
rected toward the back of the blade along the walls 
and toward the front of the blade throughout the core, 
as shown in Fig. 5. 

Another basic phenomenon, important for the flow 
through blade grids, requires discussion here—namely, 
the different behavior of boundary layers when adher- 
ing to a solid surface as contrasted with boundary layers 
detached from a solid surface, the latter usually re- 
ferred to as wakes, indicated on Fig. 6. 


Under the action of shear stresses of surface friction, 
a boundary layer attached to a surface tends to grow 
in thickness, which growth may be intensified by de- 
celeration or counteracted by acceleration under the 
effect of pressure gradients. In detached boundary 
layers, the shear stress resulting from skin friction is 
no longer active, and, therefore, the boundary layer 
thickness does not continue to grow unless subjected to 
a field of positive pressure gradients. On the other 
hand, the action of turbulence causes an exchange of 
momentum which tends to dissipate the wake. The 
rate of decay of a wake after its detachment from a solid 
surface is of intense interest in the design of turbo- 
machines, since on it depends the distance allowed be- 
tween successive blade rows. It is also a factor that 
influences the noisiness of turbomachines. In the fol- 
lowing, use will be made of the fact that detached 
boundary layers decay rapidly, while boundary layers 
adhering to a solid surface persist or continue to grow. 
The decay of wakes is a process involving relatively 
high losses of kinetic energy, the magnitude of which 
is affected favorably by acceleration of the entire fluid 
stream and distinctly adversely by retardation. Re- 
tardation of the fluid stream as a result of the decay of 
wakes, and consequent positive pressure gradients, are 
likely to occur downstream of a blade row unless com- 
pensated by changes of flow area. Typical variations 
of pressures downstream of a grid are indicated in Fig. 
6. 

A third factor of major influence upon the three- 
dimensional character of flow in moving blade grids is 
the effect of centrifugal acceleration upon the fluid 
introducing radial motion and corresponding radial 
momentum components and affecting the boundary 
layers of the moving rows in particular. As a conse- 
quence, boundary layers of backwardly sloping blades 
are displaced radially outward; of forwardly sloping 
surfaces, inward. This action is equivalent to bound- 
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ary-layer removal on some portion of the blades, re- 
sulting in accumulation on others. 


A CRITERION OF RADIAL STABILITY 


The occurrence of displacements of flow referred to 
invalidates the assumption of flow on coaxial surfaces, 
and it appears desirable to establish a criterion for its 
existence or nonexistence. The following formulation 
is offered: 

A stable flow on coaxial surfaces exists if a particle 
of fluid, when displaced radially and thereby undergoing 
a change of state of constant energy, corresponding to 
the variation of pressure of the surrounding fluid along 
the path of displacement (however, not affecting the 
meridional components of acceleration), is subjected toa 
resultant force tending to return it to its original posi- 
tion. 

Indifferent equilibrium of flow on coaxial surfaces 
exists if the resultant force is zero; instability prevails 
if the resultant force acts in the direction of further dis- 
placement. 

This criterion is applicable also to boundary layers 
along the casing and the hub.’ 

It may be expressed in differential form for compres- 
sible and incompressible flow, 


ata iis i 
Or? < or . Y / particle 


with the additional condition 0h/0r = 0 applied to the 
displacement motion of the particle. The > sign in 
Eq. (5) applies to stable condition, the = sign to in- 
different equilibrium, and the < sign to instability. 

Development of Eq. (5) leads to the following expres- 
sions for the condition of stability of flow on coaxial 
surfaces of revolution, for incompressible flow: 


0*p/dr? 2 —(3/r)(Op/dr) (6) 
or, in terms of total head h 
(Oh/Or) feta 4 0 (7) 
and for compressible flow: 
oe = + (2 i. (8) 
Or? < vyp\or r Or 


or, again in terms of total energy h 

oh 2 ) | 

— — RT| —— = (log, T) — — (keg, 3 9 
dr Fest )— g ue jee ©) 


In all instances, the upper sign applies to stability and 
the lower to instability. It is noted that stability of 
compressible flow depends not only upon the radial dis- 
tribution of total energy, as is the case for incompressible 
flow, but also upon the temperature distribution. 

Eqs. (8) and (9) are reducible to the form of Eqs. 
(6) and (7) upon introducing the condition of incompres- 
sibility. 


SPECIFIC EXAMPLES OF BOUNDARY-LAYER 
DISPLACEMENT 


Consideration will be given now to the action of 
boundary layers in the blading of turbomachines for 
the following conditions: 


(1) Fixed grids with zero circulation at outlet. 

(2) Fixed grids with circulation at outlet. 

(3) Rotating grids with zero circulation at outlet. 

(4) Rotating grids with circulation at outlet. 

(5) Unsymmetrical boundary layers upstream of the 
grid. 

(6) Flow through blade grids operating in a rota- 
tional field of flow. 


The type of boundary-layer flow occurring in the 
above grids is illustrated in Figs. 7 to 10. Fig. 7 shows 
a fixed grid with axial discharge. The boundary layer 
originally at the casing and the hub only is partly dis- 
placed to the back of the blade. Upon detachment from 
the blade a part of the boundary layer originally adher- 
ing to the wall of the casing and to the hub is inserted 
into the fluid and thereby subjected to rapid decay. 
This is a feature that commends the design for axial 
discharge from a stage. 

Fig. 8, a schematic presentation of the flow through 
a fixed grid with circulation at outlet, shows a displace- 
ment of the boundary layer from the casing to the back 
of the blade which is more accentuated than the corre- 
sponding displacement of the boundary layer of the 












BounpaAry Lavees 
AT_OUTLET.. 
BOUNDARY LAYERS 
AT INLET OF BLADE 











| ee 
aK, 
SY 
7——}—- DIRECTION OF 
Bou NDARY LAYER 
DISPLACEMENT 
SSL=“ Lecaniom oF 
Sea; ae BouNnPARY Layee AND WAKE 


OF A FIXED GRID WITH AXIAL 
DISCHARGE. 


DISPLACEMENT OF WAKE 
TOWARD THE ROOT, pe SSE 
* 


BounDary LAYeeS 
AT_NLET. 





‘ 
‘ 
\ \ ‘ : 
WN — DISPLACEMENT FLow 
«<* WN Naa - 
pe SSA cL BounDARY LAYERS _ 
‘WAKE REGION AT _OVTLET. 


BOUNDARY LAYER AND WAKE _ 
OF A BLADE Row wiTH 
CIRCULATORY MOTION AT ouTLeT 


Fic. 7 (top) Fic. 8 (bottom) 





656 


JOURNAL OF THE AERONAUTICAL SCIENCES—NOVEMBER, 





“DIRECTION OF DISPLACEMI 


——BouNDARyY LAYERS _ 
AT_OUTLET 


Bounpary LAYER awp WAKE 








WAKE REGION. 








*Bounpary LAYees 
AT INLEr 





Bouwpary Layers AND WAKE 
OF A ROTATING BLADE Row | 
Wim CIRCULATORY Motion AT OUTLET 


Fic. 10 (bottom). 











Fic. 9 (top). 


hub. A pressure gradient downstream of the grid 
further tends to displace the wake toward the hub. 

Boundary layers along backwardly sloping moving 
blades, Fig. 9, and along the rotating hub are displaced 
toward the region of the tip by centrifugal acceleration, 
as stated before, and, in the case of axial discharge from 
the blade row, an accumulation of wake fluid ‘at- 
tached”’ to the downstream casing wall may be the re- 
sult. If the fluid at outlet of the moving row has a 
circulatory motion, however (Fig. 10), the existing radial 
static pressure gradients will induce a displacement of 
the wake radially inward and, counteracting the pre- 
ceding outward motion, tend toward a more even dis- 
tribution of wake fluid. 

In discussing condition 2, (Fig. 8), it was noted that 
the secondary motion was unsymmetrical with respect 
to the plane of symmetry of the blade channel. This 


1947 


phenomenon, which is equivalent to a twisting motiog 
in the boundary layer, may originate also from causeg 
other than a field of pressure gradients—e.g., from un-) 
symmetry of the flow at inlet of the curved channel ig, 
respect to the thickness of the lateral boundary layers 
If twisting motion is imparted not only to the boundary 
layers but to the main body of the fluid as well, rela/ 
tively large amounts of kinetic energy may be involved, 
and direct and indirect losses are apt to be large. 

Considering now the action of a blade grid in a field: 
of incompressible flow in which the total energy ine 
creases with the radius, it may be shown that in thi 
case displacement of fluid normal to the coaxial stream 
surfaces is effectively opposed by stabilizing forces, 
and, consequently, secondary flows may be expected 
to be small. 

In conclusion, it should be stated that, whereas for 
single-stage units, the boundary layers generally a 
relatively thin and the effects discussed therefore min- 
ute and probably negligible, evidence from tests of 
multistage compressors, particularly, seems to indicate’ 
that the three-dimensional effects in question have prac 
tical significance. | 

The analysis of boundary-layer displacements seems] 
to offer an explanation of the observation that, par- 
ticularly for blades of short span, it is the tip section off 
the moving row rather than the root section, as might 
be assumed from two-dimensional theory, which proves 
to be critical in respect to stalling. $ 

While secondary motion introduces several factors) 
responsible for losses in the stage, the action of blade 
rows by which attached boundary layers are detached™ 
and thereby exposed to rapid decay—a process to 
which the term ‘‘metabolism’”’ might be applied with 
some justification—deserves the attention of the de+ 
signer, since it might very well be applied consciously 
to the improvement of the characteristics of multistage} 
units in particular. a 
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